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THE CAUCHY PROBLEM FOR AN ELASTIC DIELECTRIC IN A MAGNETIC FIELD 


SYLWESTER KALISKI (WARSAW) 


1. Introduction 


In Ref. [1], the problem of magnetoelasticity is treated for a body with good con- 
ductivity, for which the conduction currents dominate the displacement currents. 
It was found to be possible to reduce the resolving functions (described in the general 
case of coupled mechanical-electromagnetic fields by partial differential equations 
of the 12th order) to a simplified form, such that they are described by equations 
of the 2nd and 4th order, which facilitates in an essential manner the solution of 
this complicated problem. In addition the solutions were constructed so as to make 
use of the known equations of the classical theory of elastic and anelastic bodies. 


In Ref. [2], an analogous problem was solved in quadratures for the much simpler 
problem of a perfect conductor, for which the equations describing the resolving 
functions were also of the 2nd and 4th order. The object of the present paper is 
to obtain a solution of the Cauchy problem for the other extreme case (the first 
being that of [1]), that of a dielectric, where the conduction currents are zero. Con- 
sidering mechanical-electromagnetic fields excited by mechanical agents, simplified 
solutions are obtained in this case also, similarly to [1]. The general 12th order 
partial differential equations for the resolving functions are reduced to ordinary 
wave equations and a fourth order equation, the latter containing, in a particular 
case, a small parameter and reducing to a Volterra integro-differential equation 
with a kernel constituting a combination of solutions of wave equations. The 
equations and solutions for dielectrics are much simpler than for conductors, [1]. 


In the solution of the general problem there is an unsolved case, where the con- 
duction currents are of the order of the displacement currents. In this case, it is 
not possible to find such simple solutions as those of [1] and the present papers, 
the full solution of the problem thus remaining an open question. The idea of the 
construction of the resolving functions is the same in the present paper as in Ref. 
[1]. The notations are also so chosen that the system of equations is apparently 
identical with those solved in [1]. The similarity is apparent because the same sym- 
bols correspond to different operators and parameters. In spite of the identity 
of the methods for constructing the resolving functions with those used in [1], the 
basic argument will be repeated here in view of the two papers being written in 
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parallel, and because that argument is of essential meaning for the principal stage 
of the solution procedure, that is, the construction of the simplified equations 
which are of a different and much simpler form than those of [1]. 

In every case where the conclusions of [1] are preserved we use them without 
proof (for instance, in the case of the construction of the operational resolving 
functions for anelastic bodies). 


2. The Resolving Functions 


The linearized equations of a mechanical and electromagnetic field for an iso- 
tropic real conductor in an initially uniform magnetic field H have in the general 
case (that is where the conduction and displacement currents are of the same order) 
the form: 


( a AT ease OK eu—1 dO [du 
rot h eet Fey: ce alan): 
rot pee OL 
c Ot 
ap Pea Se 
(2.1) Alara eine or \} 
Cu ie 
ose = GV2u-- (A+G) grad divu-+—— GxH)+ 


1 Sct Lt 0 [ou \ 
lamas D (Fe) + gta ) ar (ap), 


and 

(2.2) divi .02" “dive Dp... 

where 

(2.3) pa fete al, 
ce Ot 


On the basis of (2.1), (2.2) and (2.3) we have at the same time 
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(2.4) div j= 


In the case of dielectrics these equations are simplified. We have, [3]: 


(2.5) j= OP peer 
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Then, the system of equations takes the form: 


e OF  eu—l oO /du 
MS. c Ot ce artar<H): 
rot E= — fae 
C Gt 
(2.6) 
O7u 5 : 1 OE 
oa GV2u-+(A+G) grad div u+ yee (ice (zr) <H)+ 
Ub nee ren 
[ aa ae (Sapa, 
and 
(2.7) dv h—=O; dizD=0; 
Hence 
: sae 1es/,0u 

(2.8) div E= = div | ar <n| ; 


The introduction of the general form of the resolving functions in the system 
of equations (2.6) would, as was shown for an analogous system of equations in 
[1], enable the reduction of this system to separate partial differential equations 
of the 12th order for the resolving functions. However, proceeding as in [1], and 
introducing certain simplifications consisting in rejecting very weak couplings, 
the problem can be simplified considerably and reduced to the solution of equations 
of much lower orders, for the resolving functions. 

With no limitation of the generality of the solution we can assume H = H; = H. 
Then, eliminating h from (2.6), we obtain 


2 <2 3 
rotrotE}+"° gen eae lee xH) =0, 


Ot G or? 
OE. 
Kee Gs divu = AA i ty Se 
j Q OX, at 
(2.9) 
1 re) OE, 


[ Kity——-(A4+G) 30 divu = oP 
where 
f a Fa|'+ ee =| av’, 
0? 2072 
(2.10) K= 7a HV’) 
pe—l G c : 
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Making use of the relation 


(2.11) rot rot E = graddivE—V’E, 
the system of equations (2.9) can be rewritten [bearing in mind (2.8)] in the form: 
( eu—1 Ou, ,, me-l_ @ (Ou, Om\He 
| pes ec 08 a! I~ Ox,0L 0X, age |e 
eu—1 Pu, pe—1l CP OU, 6a; Vile © 
(2.12) 4 NED ec of eu = OXg0t\ OX, Ox, ) € e. 
| eu—1 ole Cu, . Ou, \-He 
a =): 
sie eu OX, 0t \ Oxy 0X2] € : 
( AG 0 ~ cay Bod OE, 
| Kyu,— et divu Se ae 
A+G 0 OE, 
(2:33) 4 Ri are Ox, divu = an Py—y a? 
| A+G 0 
[ ae aa ax, divu = —Ps3, 
where 
2.14 a 272 
(2.14) N= a —ooV . 


In the operator K, (2.10) we have H*u(ew—1)/(4ec”) < 1 and 
be disregarded in relation to unity. Then 


(2.15) Roe 


Eliminating £,, E, from the last three of the equations (2.13), and 
the first three, we obtain: 


this term can 


substituting in 


( NKu,— 6 i Ndivu— a viT a 
ara 4 ce et ee 
(2.16) Ngee w Ndivu— eo zs = L 
PMN gf BeBe) HE Le, 
[ Kuz are. “ divu =—Ps = a 
ca Nie aaj tq ee 
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Thus, similarly to [1], the original system of equations, for which the resolving 
functions were described by equations of the 12th order, has been reduced to the 
system of equations (2.16), which, after disjoining, reduces to an equation of the 
10th order and the wave equation (2.17), from which E, is determined after deter- 
mining u from the system of Eqs. (2.16). For the excitation force field, it is as- 
sumed that these forces admit the operation N. The system of equations (2.16) may 
in what follows be reduced to equations of order lower than the tenth. To this 
end, let us introduce, similarly to [1], the following set of resolving functions: 


Ov, OF, | _ Wi, oP, OW, 


ae eine Cx Ok. cia aiar ue Fe mae sem Se 


Substituting (2.18) in the system of equations (2.16) we obtain: 
7) AG oo oN, A+G OF 
( ee eames N a = (Rr 5; yar Vi|Pe= Pi, 


OX, 
é A+G ie Cee 
ee Sree Mss ar Vi | was, 
(2.19) 1 x, (8 ; vi] 2a 5 ee alers dl Ayes 2) 2=P; 
0 Ge A+G @ =e 
ea 5 iV 44 (« meee :)¥|=P2, 
where 
4 oO? akg 
Ch 5 
Se eee axe ~~ Ox2” 
pe ev Se (gu—1)? . 
Bay iw Seeoe re ~ 4gcec?o te 
Zi (ease EAE) ers 
ae gan Oe Ao eu 
Representing P; in the form 
OXq- 0X3 OX BOX, OX 
be ee ae * fae Oe ea . P. a ; 
Ce) oes Ox, OXn. P, OXe us Ox, Bet Ox, 


where X; is calculated, after transformations similar to [1], according to the equa- 


tions 
Decent 0 I 
(ez if i} (? on thf go In Jabs, 
| ‘ Gomes 
0 1 
= ; Ps dé,dé,, 
(2.22) 4 Xs [fe Ox, ln r i ae a] Ase 


7? = (x,—€,)*+(%2—¢:)"; 
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the system of equations (2.19) can be rewritten thus: 


a 
(2.23) (+7-57¥3) Yo = Xo, 
( band la AG as ie 
(2.24) 
Wer G. AGO 
esa ame ee cs YY, a Xo. 


The initial conditions for the system of equations (2.13) will be assumed in the 
form: 
Ou, (1, Xo, X3; 0) 


Ot == £\%4,.¥25 Xa) 


f U(X, X25 %X3, 0) =f, X25 X3)5 


(2.25) 4 
{ E,(%1, Xo, X3, 0) i 
where f,, g, are assumed to be of the C? class. 

The functions u, in the system of equations (2.16) satisfy the following initial 
conditions: 


OE, (%1, X25 X3, 0) 
Ot 


mY 


Oth; Cees Xe 0) 
— — = 8, (%1, X25 Xs), 


[ Uy, (ae X25 X 35 0) = fi,(%1, Xo5%s,)> 


OPuj(X1, X2, X3, 0) = A+G 0 


div f-+-a? BV + oP, (x1, X25 X3, 0) 


(2.26) ¢ ot? @ Ox; (i= 1,2), 
| -_ 
OPuj(X1, X25 Xs, OEP ; AG ¥e 2 1 OP, X2, Xs, 9) 
| or? OaeOr; div g+aV'a-1 Q Ot 
t @= 1,2). 


The functions yw; satisfy, according to the equations (2.21) the initial conditions 
( P (1, X25 X%3,0) = bj (%1, X2, X3), 


0 ’ 2 0 2.4 
PO = X3 Dis mae Doe, Xs, =) (k = 2) 
(2.27) o* ( > > 0 Xi 
{ ¥ Jj =e X3 ) — by (%4, Xo5 a), 
8Us. 1 = 
0 suet X3,0) (Ka Xan Xa) (j= 122); 


where b,, d, are determined from the equations] 


0 1 0 1 
(2.28) b= [{(% ax, In ; +f, ax, In Jats, 
i Qy 
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0 1 0 1 
b, = 
a { [ (4a - r fray, aE r Jab 
Q, 


by = i fad, 
0 1 0 1 
Cr | 
1 {{(« Ox, In r i 82 OX In r Jie, 
25 
ro) 1 0 1 
d, = 3 ein ot 
2 {i fl (2 OX 2 ms! OX a r Jaede, 
2, 


ds = | sodés. 


(2.28) 


The functions 5;, d,; may be calculated by means of the functions h,, d, on the 
basis of the relation (2.26). 

Let us now disjoin the equations (2.24). We shall introduce new disjoining 
functions ®,, ®,: 


( ALG @ AG 3 
ge bee Ss mee 
(2.29) 
A 
[as vio. (R Nv). 


Substituting (2.29) in (2.24), we obtain: 


(2.30) [Gaerne se A+G @ “| (AG)? 
e Ox? 0? 


NVI, |e ei mi EGS 


These equations are apparently identical with those obtained in [1], it should 
however, be borne in mind that the symbols, which are-identical in both equations, 
denote different expressions. The Eqs. (2.30) are simpler, in principle, than the 
corresponding equations of Ref. [1], because the difference between the operators 
K and kK, vanishes. 

The functions ®; satisfy the initial conditions: 


(2,6, % %)— ete g = 13), 
OD Cia XsiXe0) OD as eon Xe) 
iy Tete = bit ty xy),  — Sage? <0, 
OD ks, X50) 07D, (X41, Xo, X3,0) 
* ae ° = d,(X1, X2, X3), ° a o =0, 
AO,(%1,; X25 X3,0) 


(2 31) att 7 bya, X25 x3) — 
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(2.31) r ane oe SEMEN ok +(av+ Ze gag) Xay Xa) 
Dee TOR 
Meee: ee Ay 1X1, Xa, %3)-+ 
+ (aves ae Sag} Xe, ye Feu Xa» Xs)> 


OD, (Xs Xo, X3,0) 
dt® 


The aim is thus reached. Similarly to the Ref. [1] for conductors, in the case 
of dielectrics the general equations of the 12th order for the resolving functions 
reduce to equations of the 2nd, 4th and 6th order ((2.17), (2.23) and (2.30), respec- 
tively). This essentially facilitates the solution of the problem, but the solution of 
the 6th order still remains difficult. 

Similarly to [1], the equations reduce in the plane problem to the 4th order. 
Then u, = E; = 0, and we obiain, on the basis of (2.23) and (2.24): 


= d3(%1, Xo, Xg)- 


(2.32) [Riv Savi)? a (z- owva a OF 
4 

To avoid the solution of 6th order equations we shall consider the possibilities 
of simplification of the equations obtained for the resolving functions by rejecting 
certain small physical couplings of higher orders. This will enable the problem 
to be reduced to the solution of wave equations and one equation of the 4th order. 
The simplified equations thus obtained for coupled mechanical-electromagnetic 
fields in dielectrics are much simpler than analogous simplified equations of 
Ref. {f]. 

In the Ref. [1] two variants of simplification were considered. In the present 
paper only one variant will be considered—the other variant of the Ref. [1], 
concerning the dissipation. 


3. Approximate Equations 


The quantity v7 in the equations (2.23) and (2.30) can be rejected as a small 
quantity of a higher order. Similarly to [1], we shall show by means of an anal- 
ogous one-dimensional example that »y is of the order a3/c?, and may be disre- 
garded in relation to unity independent of the time ¢ at which the process of wave 
propagation is examined. Putting »7 = 0 in the Eqs. (2.23) and (2.30), we shall 
obtain the following simplified forms of these equations: 


02 
G.1) [VK-y aN) Bs 
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A+G A+G @ QA+G) 03 
ays Nii K— 2\/ e_ =. Fie Jets < 
Be { e vil(k g oa) ee vila, Ae SA3)- 


The equation (3.2) can also be rewritten thus 


(3.3) NO102®; = Xx; (i a 1,3); 
where 
ig 2072 0? 
Cy Oz ayV’, Ll, = ap RV" = K, 
(3.4) 
| eG a ié 
ay — 0 as —= <= 


Thus, the equations of the problem have been reduced to the solution of wave 
equations and one equation of the 4th order instead of the equations of the 2nd, 
4th and 6th orders, in the case of 7 40. 

The Eqs. (3.3) are of the 6th order, but they are expressed by means of a triple 
product of wave operators which can be integrated in a simple manner by solving 
three wave equations succesively. Denoting 


(3.5) 0,022; = 9;, 


we obtain from (3.3) the equation: 

(3.6) NO; = X;, 

where the initial conditions follow directly from (2.31). 
After determining 9; from (3.6) and denoting 

(3.7) 0.8; =T;, 

L; 


(3.8) OL; = 9; 


can be calculated from the equation 


with appropriate initial conditions for I’,. 

I, being thus determined we calculate, from the wave equation (3.7) the functions 
®, sought-for. Thus, finally, the problem is reduced to the known solution of 
the wave equation (2.17) and the wave equations (3.6), (3.7), (3.8) and an equation 
of the 4th order for yw, (3.1). 

The latter can be written in the form of a Volterra integro-differential equation: 


G9) =f ff ees a] HG 0-1 GaVEvAG 9 [dar 
oa 

where G is Green’s function for the equation 

(3.10) OoLeG = o(x—§, 4), 

where 


Bese 


246 S. Kaliski 


The Eq. (3.10) will be solved, as before, by solving successively in a known 
manner two wave equations. The formulae which are known will not be cited 
here. ; 

With the assumptions on X; previously made the Eq. (3.9) has a convergent 
solution. If the original field H is small, y may be treated as a small parameter 
and we can confine ourselves to the first approximation only. It will be shown 
now, on a simple, one dimensional case that vj can be made equal to zero in 
the system of equations (2.13), which is equivalent to the rejection of a quantity 
of the order a?/c? in relation to 1. 

Similarly to [1] let us consider the case where u, = u, = E; = E,=0, Py = 


= P, = P,; = 0. The initial conditions for E, will be assumed to be homogeneous 
and those for uv, to have the form 


OUy(x,0) 


(3.11) Us(x.0)-—= @&), ar 


0. 


Then the system of equations (2.13) takes the following form, assuming that 
the solutions depend only on ¢ and x, = x 


2 8) Om uy 
7 fe 2 =) Gx OR” 
(3.12) a2 ess 
Oe ae) ot ° 
where 
(ye ae ae ee pe rey 
be amore ; eo 4aco 


where » may also be considered to bea small parameter; therefore the above 
equation can be solved by means of the iteration method™, by determining u, from 
the second equation, then £, from the first, u, from the second ete. 

To make an appraisal of the right-hand sides of the first of the Eqs. (3.12), 
it suffices to substitute u, of the zero approximation and to set equal the two 
members. Hence their ratio is obtained. 

The zero solution for u,(x,t) has the form: 


G.14) ua(xt) = 5 [p(x—aat) + peta). 


Substituting (3.13) in the right-hand member of the first of the equations (3.12) 
we have: 


Ou Ou 1 mt a 
(3.15) De U ae = >, [(—yay-+nay)y (x—agt) + (ya,—na3)p (x+-apt)]. 


“) The strong solution may also be found without difficulty but it is not necessary. 
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Setting equal the coefficients of y’’’(x—a,t) and g’’(x-+-a,t) we have 


(3.16) ya,>743, 
because, making use of (3.19), 

(3.17) NY Ae leer Ny 

pe EC 

Hence 
: C2 ce 

3.18 Zee 

(3.18) Bue ae >I, 


which confirms the assumption. 


4. Resolving Functions for Anelastic Bodies 


The considerations of the Sec. 2 are valid formally also in the case of anelastic 
bodies. The solutions being identical with those in the case of Ref. [1], we shall 
give final results for the resolving functions. As an example, a Boltzmann 
solid will be considered. It is assumed for simplicity that the initial conditions 
are homogeneous. Subjecting the system of equations (2.13), with additional re- 
laxation terms for a Boltzmann body to the Laplace transformation and intro- 
ducing, as in Sec. 2 the transforms of the resolving functions, which will be called, 
for brevity, operational resolving functions, we obtain: 


(4.1) (R+yp’V?) = X), 
Gee i ALGO? GsEG) a. 20 / ae 
a, a we op a = 13 
oe (R 0 yi) (« 0 =| e ae Ox} it 
or, assuming 7 = 0, 
(4.3) (NK+-yp°V2)P2 =aXoe 
(4.4) NOOO, = Xx; (i = 1,3), 


where the barred operators denote the operators defined in the foregoing section 
except that 0/dt is replaced by p, A, G— by A, G, and 


(4.5) Ai =R@).9 1G —G—O), 
where R(p), Q(p) are the transforms of the relaxation functions appearing in the 


equation 


(4.6) 4 = Adivud,+2Gen— ‘| [R(t—1) div u(t)d,+20(t—1)e4(a) dt. 
0 


The solvability of these equations is, depending on the relaxation functions, similar 
to the theory of anelastic bodies. This problem will not be dealt with here. 
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5. Conclusion 


The present paper and Ref. [1], enable the construction of a solution for com- 
bined mechanical-electromagnetic fields in a relatively simple way for real conduc- 
tors and dielectrics. 

In these cases, the general equations of the 12th order for the resolving functions 
can be reduced to equations of the 2nd and 4th order, for which the solutions 
of the Cauchy problem are in most cases known in the theory of elasticity. This 
considerably simplifies the solution. In Ref. [2], similar solutions (the equations 
of the 2nd and 4th order being solved in quadratures) are obtained for a perfect 
conductor. 

If the conduction and displacement currents are of the same order the above 
solution cannot be directly used, and that problem requires separate treatment. 
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Steric Siz Gl Zremame 


PROBLEM CAUCHY’EGO DLA SPREZYSTYCH DIELEKTRYKOW 
W POLU MAGNETYCZNYM 


W pracy rozwazony zostal problem Cauchy’ego dla odksztalcalnego dielektryka 
sprezystego i niesprezystego znajdujacego sige w pierwotnym, statym polu mag- 
netycznym. Podobnie jak w [1] udalo sie w przypadku wymuszen polami mecha- 
nicznymi rozbi¢é ogdlne réwnania czastkowe 12 rzedu dla funkeji rozwiazujacych 
na szereg r6wnan falowych i jedno rédwnanie 4 rzedu, ktére z kolei sprowadzono 
do réwnania rdézniczkowo-calkowego typu Volterry z jadrem bedacym kombinacja 
rozwiqzan rownan falowych. Tak wiec zagadnienie zostato sprowadzone w zasadzie, 
do kombinacji rozwiazan znanych. Rozwiazania otrzymane dla dielektrykéw 
opieraja si¢ ze wzgledu na podobny charakter rownan wyjsciowych na analogicznej 
konstrukeji funkcji rozwiazujacych jak w [1] z tym, ze rozwiazania dla dielektry- 
kow ksztaltuja sie znacznie prosciej anizeli dla przewodnikéw ze wzgledu na inny 
charakter uproszczen rownan wyjsciowych. 


Pesrme 


SATAUA KOWMM DJIA YIPYTAX DTAMSJIEKTPHKOB B MATHUATHOM IIOJIE 


Paccmarpusactca sayjaua Koum aia gedopmupyemoro yupyroro mu Heyupyroro 
AMSNEKTPUKA, HAXOTAWIETOCA B IICPBKUYHOM TIOCTOAHHOM MarnuTHOM rome. Ilo- 
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qo6HO, Kak B padore [1], ymanocb, B cilyuae BHInyUeHuii MexaHMueCcKMM OJIem, 
pa30uTb oOnjee ypaBHeHHe B UaCTHBIX MpoUsBOAHEIX 12-ro nopamKa WIA paspema- 
FOINMX (PYHKUUH Ha PA], BOUHOBBIX ypaBHeHHH wu OWHO ypaBHeHue 4-ro mopayKa, 
KOTOPO€ B CBOIO OUePeb CBECHO K HHTerpo-zuNpdepeHuMabHOMy YpaBHeHHIO TMA 
BossTepyi, C AUPOM ABIIAIONWIVMCA KOMONHAaUMCH peuieHHii BONHOBBIX ypaBHennit. 
Takum oOpa3om 3ayjaua CBOXMTCA B OCHOBHOM K KOMOHH@IIMM H3SBECTHbIX pelleHnuit. 
PelicHHA MONYUCHHBIC WIA DMIIEKTPVKOB, BBUTY moROOHOrO pemleHuA MCxOJ- 
HbIX YpaBHeHHi, OCHOBLIBaIOTCH Ha aHaJIOrMUHOM KOHCTpyKUMM paspemlarompx 
dypKuni Kak B [1] c TEM, UYTO PCUICHHA NA WusIeKTpUKOB MO cpaBHeHHIO 
C NpOBOAHHKaMH, MUMeIOT OoNee MpocTy!lo MopmMy, BBULAY MEHOTO xXapakTepa 
YNpPOWeHHH UCXOQHBIX ypaBHeHHi. 
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PROCEEDINGS 

OF VIBRATION PROBLEMS 
WARSAW 1961 

VOL. 2, No 3 (8) 


THE VIBRATION OF A CYLINDRICAL SHELL OF FINITE LENGTH 
WITH A SUPERSONIC INSIDE FLOW 


JERZY NIESYTTO AND ZBIGNIEW SEP (WARSAW) 


Introduction 


A number of rocket engineering problems—those of rocket and ramjet nozzles 
for instance — consist in determining the conditions of self-excited vibrations 
(flutter). The physical model for these problems is a cylindrical or axially sym- 
metric shell of finite length acted on by a gas flowing inside the shell. The way 
of fixing the shell (the boundary conditions) and the so-called structural damping 
is determined by the design features and the so-called material damping — by the 
kind of material used. Each one of these factors reduces, in general, the possibility 
of vibrations, but the appearance of them is not entirely prevented, however. 

In the literature there are some works concerning the flutter problem of infi- 
nitely long shells flown past a supersonic flow outside the shell, [1], [2], [3], [4], 
[5], [6]. Material and structural damping and also the anisotropy of the material 
are taken into account in [5] and [6] alone. For the analysis of the flutter problem 
it is important to have the high degree of accuracy of determining the gas pressure 
on the shell, provided that the possibility of actual computation is not impaired 
by mathematical complication. 

In the recent works [7] and [8], considerable simplifications are made for the 
aerodynamic forces, or only a special vibration is considered (axially symmetric 
vibrations). 

In the present paper a flow inside a shell is considered on the basis of the 
linearized theory of flow and the equations of the shell are deduced according to 
the general bending theory of shells (cf. [9], [10]). The flutter problem of a shell 
of finite length requires the satisfaction of the boundary conditions for the potential 
of flow and the equations of motion of the shell. This complicates the problem 
from the qualitative and numerical point of view. 

The solution presented enables actual numerical computation. The computation 
algorithm proposed may be adapted for computation by means of electronic com- 


puters. 
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The method of analysis of the influence of each particular structural and aero- 
dynamic parameters (diameter, thickness and length of the shell), the material dam- 
ping coefficient, the flow velocity of the gas etc. — is based on the Neumark-Nyquist 
method which enables convenient graphical representation of the results [13], 
[23]. 


1. The Statement of the Problem 


Let us consider a cylindrical shell of finite length the radius of the middle surface 
being R and the wall-thickness A (cf. Fig. 1). 


Fig. 1 


The flutter phenomenon of the shell will be investigated for a gas flow through 
the shell of which the velocity of the unperturbed supersonic flow is U. The ma- 
terial damping of the shell will be taken into consideration by assuming the Voigt 
model, similarly to the Ref. [5]. The stress-strain relation can be assumed in the 
form 


0 0 
CC (a+ iver} A +2(urbee 5) Ex» 


) ) 
es (at ar|4 (a Me sr) 


0 
Tx@ a (ary Exe py 
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where A= ¢,+6,, €,, &—unit elongations in the directions of the axes u, v 
(Fig. 1), €,,—change of the angle between the axes u and y, 
Ey E 


ema y eS Che 


where: A,—volume damping coefficient, ~,—distortion damping coefficient, E— 
Young’s modulus, »—Poisson’s ratio. 

On the basis of the general bending theory of shells, [9], we obtain the equilib- 
rium equations in the form 


( ou ow 
(24-20) — theese + Ot ae $1 


O'w O7u 
—c(A Pryeaa 3 f OP aah — RO. aa = 0, 


9 
(1.1) (axaacaaz Habe +04 3H SE + At get O42) 5™ 


9 08 w 0?y 
a ie Le ia 209 Ro, A 0, 


Ou 
(e.- : cesap ale a uer-s Gee 


are Paes, 


Otw o*w 
— 20 


ee o4w 
TATA ayaa sags * ait 


oie cee o?w ou , ov 
5¢ A 5a 2c?(A- Wace fae agi 
oC? R* 
+(40°)(1+2p)wt Rosa — sp = 00) 
where 
2 he 
Coe i= T2R2’ q+, | 
u,v,w displacements of the middle surface in the direction of the respective axes 
(Fig. 1), 


p difference of gas pressure of perturbated flow and unperturbated one 
inside the shell, 

0, mass density of the shell, 

o@ mass density of the gas. 


() We assume that the pressure difference between the unperturbated flow internal and 
the external (immobile) gas can be disregarded as well as the influence of the radiation on 


damping. 


2* 
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In his theory of cylindrical shells Vlasov, [9], rejects the terms containing c* and 
c2. Fliigge, [10], rejects only those containing c’, the terms containing c? are 
retained. In the present considerations the terms containing c* will be the only 
ones rejected(’). 

We shall consider a shell rigidly fixed at one end, the other end being free (Fig. 1) 
The fixed end conditions are therefore such 


0 
Gea) Ly == = - 0 for. 3x0 
Ox 
For the free end we have: 
La= 0) 2 0 for x=I, 
(1,3) Smey 
NM Os tm a =0 for x= 1 


Fig. 2 


The conditions (1.3) may be expressed in the displacements u,y,w: 


Ou ov ne w 
2 = 
(A+2p) eo ae AW er(A-- 2 ieee = 0, 
Ou 0?w 
2 
lb a r(1-+-3¢ ee ax = ey 0, 
(1.4) ia p 
(A+ 2) ie (A iy i =0, 
dy 
0*u Ou On 4 ow 
| ee 
[ (A+2u) A323 Haga (A4 axe —A+2 ee ery Koper y 


The inside gas flow under consideration is acted on by the vibrating shell. The 
equation of the potential ® of a linearized flow has, in cylindrical coordinates, 
the form: 

OD OD oD 


eS 2AD = | 72 
Oe a ot? ape oo Oxee 


(*) The influence of the c and c, is appraised numerically to some extent in Ref. [14]. 
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where 

o® 1 o@ o® 

beg 71 ON Gh RE ok ey oe 

Ox o( x )s r Op OV os or ee 
a — velocity of sound in gas, U — undisturbed flow velocity, V,, V,, V, — velocity 
components in the respective directions. 


The potential ®(x, r, g, t) will be sought for in the form 
DCT. ¢, 1) = K (YO, o, 2). 


From the boundary condition of flow on the wall of the shell 


Ow . Ow e = 
1.6 = 3 


relating the potential with the displacements and assuming 


eee 


after some simple manipulations, we obtain 


1.7) G19.) = Free Re O[ G HU GI: 


Hence, assuming e to be constant, we obtain 
pia RQ Cow 02w ,0°W 
OE en cre: FECA ae nae Bota ae 


If in the boundary condition for the potential the 
fact is accounted for that the perturbation potential 
of supersonic flow, M, = U/c > 1, we have on the conical surface [x = x»>—(X%)/R)r] 
with the apex at the point x,, that the perturbation potential has the form (Fig. 3) 


(1.9) P=0. 


Fig. 3 


2. The Solution Method 


Particular solutions of the system of equations (1.1) will be sought for in the 


form: 
u = Aellaxtnp+et) , 


(2.1) y = Beilaxtnp+or) | 
w= Ceilax+npt ot) | 
On the basis of the Eq. (1.7), we shall search the potential ® of the form: 
(2.2) Dear, Q; j= K(r)eiextne+ en) , 


} ow ow 
(*) In further considerations we assume V, = —y7 + 5, U 
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Substituting (2.2) in the equation of the potential (1.5) we find that the function 
K(r) satisfies the Bessel equation 


pe 
(2.3) KO RO -(ne TKO) = 
where 
2 2 QU 
(2.4) m = a1 a lise a = 


The general integral of the Eq. (2.3) is: 
K(r) = ¢,N,(imr)+col,(imr). 


Hence, in view of the boundedness of the potential ® for r=0O we obtain 
Cy — 0, or 
(>) K(r) = c2l,,(imr), 


where I, is a Bessel function of the mth order of an complex variable. 

On the basis of (2.5), we can determine the ratio K (R)/K’(R) appearing in the 
relation (1.8). 

Substituting (1.8) and (2.1) in the system of equations (1.1), we obtain a system 
of homogeneous linear equations for A, B, C,: 


(2.6) a;14+4;.B+a,,C = 0 G=31,2:5); 
where 4,1, @;2, 43 are 
Ay, = 07(6,+ 6,01) +n? +c?) (uy +p.wi)— R%0,0%, 
Ag = an[(Ay+ Ya) + (Ar+ e207), 
Ay3 = —la[A,+A,wi+ c®a?(0,+0,0i1)—c?n? (uy + uwi)], 
Ag. = an[(Ay+ fy) + (Ag+Me)oi], 
zo = (1+3c?) (uy+Mwi)a?+n?(0,+0,wi)— R’o,0?, 
(2.7) { az,3 = —in{(O,+O,0i)+07c? [(Ay + 3uy)+ (Ae +3us)oi)}, 
a3 = ian®c?(u4+paw@i)—c?a°i(6, +0,wi)—ai(A,+ A,wi), 
3,2 = —in{c?a?[(A,+-3uy) + (Ap +3uug)@i]+-6,+6,wi} , 
Ag, = —-c?(0,+0,wi) (a? +n") —2c?n?(6, +0.wi)+ 
K(R) 


2 TN ey BY 2. p2 
+ (1+c?)(0,+6,wi)— R*e,w?— R iK'(R) 


0(w*?+2Uaw-+ Ua?) = 0, 
where 
0, = Ay+2,, 92 = Agt+2ug. 


In order that the system of equations (2.6) should have a non-trivial solution 
the characteristic determinant W of this system must be zero or, in other words: 


(2.8) W = Det (a, ) =0. 
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It is easy to observe from (2.7) that a; =a4;; (j= 1,2, 3). The determinant 
(2.8) depends on a and w; therefore the Eq. (2.8) has the form: 


(2.9) W (a,w) = 0 (°). 


The Eq. (2.9) contains two quantities of interest. These are: a — the deformation 
wave length and w — the angular frequency of vibration. To determine them the 
boundary conditions of the shell must be used (cf. § 3 B). 

It is impossible, in general, to find directly the roots a; of the Eq. (2.9) in terms 
of w,U. The establishment of two equations each containing a; or w only will 
be the object of a further part of the paper (Sec. 3B). 


3. The Boundary Conditions 


A. The Boundary Conditions of the Potential. The potential of flow ® [cf. (1.8)] 
must satisfy the Eq. (1.5) and the boundary condition (1.9). 

It can easily be shown that if ® = O(a,) satisfies the Eq. (1.5) and (1.6) for 
r = R, the expression 
will also be a solution. The constants Dj; are choosen to satisfy the condi- 
tions (1.9). 

Substituting for this purpose (3.1) in the boundary condition (1.9) we obtain, 
on the basis of (2.2) and (2.5), the expression 


(3.2) [®] = (ee ) = = 3 Dik (im, rei; (x)—xor/R) — (0, 


Expanding the functions J,(im,r) and e#j*4-'®) in power series of r we collect 
like powers(®) of r. Setting the coefficients equal to zero we obtain an infinite 
linear system of algebraic equations in D, 


f DAG) =U; 


Yo, [a,(a,)b.(a;)+a,(a,)b,(a,)] = 0, 
(3.3) i 


D>, D, [ay (4, )Bm(G;) +42 (0))By a (Gj) +--+ +n (4;)bi(a,)] = 0, 
be 
(*) It can easily be shown on the basis of what is known as Picard’s minor theorem 


(cf[11]) that the equation has an infinite number of roots ay (j = 102.) 
(®) The process of grouping these terms is admissible by virtue of the Mertens theorem (cf. 


(cf. [12)). 
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where a,(a,;) and b,(a;) are coefficients of terms of expansion in power series 
of the functions /,(im,r) and ei@j%(—7/R), 


The expressions “ 
“= ay! Achaea Tom), 
v=1 
(3.4) t= Dy Beilayx+np+ ot), 


v=1 


w= oe Cerri tet), 


v=1 


and the relation (3.1) satisfy the system of equations (1.1)(®). The constants A,, 
B,, C, (which have not yet been determined) satisfy the system of equations (2.6) 
The constants A,, B,, C, will be determined by means of the boundary condi- 


tions of the shell. 
B. The Boundary Conditions of the Shell. On the basis of the first two of the 


Eqs. (2.6), A, and B, can be expressed in terms of C,: 
A, = f@)C,, B, me g(a,)C,. 


Hence, by virtue of (3.4), we obtain: 


ae Dy fa,)C,eterrt ne + er) , 
v=] 

(3.5) 4 i? == » g(a,)C,eia,x+ne + ot) 
v=1 


| w= Dy Cela tae ten), 
v=1 
Substituting (3.1) in the boundary conditions (1.2) and (1.4), we obtain eight 

equations 

Pomona 60 

(3.6) on sou 
d», (4)C,=0,  d/a,C, = 0; 
v1 v=1 


f » [(0, +6,ci)a, fla,) + (A+ A,ci)nig(a,) + 
+(A,-+ Aoi) +c?(6,+6,wi)a} Ce! = 0, 


» [(O,+ O.wi)as+ (Art Agi )n? + (0, +-0,wi)a,if(a,)+ 
(3.7) 4 +(A,+ A.wi)nig(a,)]C,e%! = 0, 


(®) See also Sec. 5. 
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» [f(a,)ni+¢(a,)(1+3c*)a,i+3c?a,n]C,ei%! — 


2 (r+ O.0i)a5f (4,)— (Uy + Main? fa,) + (A+ Agwi—py —fe0i)na,g (a,) — 
— (9,+6,ci)a,i— (A, + A,wi)a,n?i]C,e4! = 
Making use of (1.6) the constants D; may be expressed in function of C;: 


_ 7 (@+aU)o 
(3.8) D,= So a 7(im,R) 
We introduce the notation C; = f!C;. 
From (3.6), (3.7), (3.3), (3.8) we obtain an infinite system of equations in the 
constants ere 


rare? 53 aa 

al ~ © a, ~ 

Qi GBE =O, DG =0, 

Dj (O.+9,01)a, f(a,) + (Aa+ Aaeri)nig(a,) + (A+ Ari) 
v=1 


+e(0,+0,oiab] 5 Cet te=\Q) 


> [O:+0,0i)a8-+ (Ay+ Aacoin®+ O, 4 Ogevi)a,,if(a,) + 
+(A,+A,wi) nig (a) ei! ='0; 


bse [f(a,)ni+g¢(a,)1+3c*)a,i+3c%a ny, el = 0, 


v=1 


Oy [(0,+0.wi)a5 f(a,) — (uy + o2@i)n? f (a,) + (Ay t Agwi— py — pani )na,g(a,)— 
v=1 


1 
— (6,+6,wi)o3i— (A,+ A,.wi)a,n® ily Cel! 03 


= m,1,(m,iR)v! 
(+ Ua,)o 


“mT GontRypt WC bA(4) +an(ay)br()IC, = 0, 


v=1 


SPE laa ba(a)+44(a br a(a) +a bale = 0, 


v=] 


ee nae —0, 
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Let us denote the coefficients of és in (3.9) by d,,, where & is the number 
of the row and uw — that of the column. Let us write the principal determinant 
of the first n equations of n variables (3.9): 

dy dy» dis eee di, 


dy 4 do» ee 


(3.10) 2"! — Det(d;,,). 


dy dy 9 dns eee dn 


It can be shown that the determinant (3.10) is normal, that is there is a limit 


(3.11) lim Det (d;,,). 


n->co 


The wavelengths a, and the angular frequencies in which we are interested will 
be determined from the relations (2.9) and (3.10) 


(3.12) W (a,w) = 0, 
(3.13) Detig )==0, 


Expanding the functions W(a,w), Det (d;,,) in power series in a, we obtain 


(ine — 
(a,@) = > Soya”, 
(3.14) | es, 


co 


Det (d;,,) = S°5,(o)", 


v=1 


or, taking a finite number of terms (,, 6, in the expressions (3.14), we have 


By 
W (a,w) + W,(a,@) = soya” == Q, 
(3.15) cae 
Det (d;,,)~ Wo(a;,00) = YS (@sea,--.M49, 1), =0 (j=1,2....,6). 
v=1 


Using the Fermat or Silvester method, [18], the quantities a; will be eliminated 
from the system of equations (3.15), thus obtaining a system of equations where 
each equation contains either a; or w: 

Bs 
(3.16) Dd 4,U)o" = ; 


»=0 


: Bs . 
(3.17) SE, (U)a} = 0. 
v=0 
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For an analysis of the influence of each particular structural parameter R,/,/, 
0,0s,A,,44, etc. on the possibility of appearance of flutter the Routh, Hurwitz or 
Michailov method can be used [21], [20], [19]. These methods, which consist in 
determining the sign of Rew [in our case we have Imw >0 in view of (2.1)] 
would require tedious computations. Therefore, to determine the region of ap- 
pearance of flutter depending on mechanical and aerodynamic parameters the 
generalized Neumark-Nyquist method will be used (See Sec. 4). 

It should be observed that the application of the Hurwitz criterion to the Eq. 
(3.16), which has complex coefficients E,(U), requires the Eq. (3.16) to be repla- 
ced with an equation of the order 28, with real coefficients. This equation is 
obtained on the basis of (3.16) by multiplying it by an equation in which the 
coefficients E,(U) are replaced with the conjugate coefficients E,(U) (cf. [23]). 
In the latter case; in order to avoid cumbersome computation, the Schur criterion 
for polynomials with complex coefficients, [24], analogous to that of Hurwitz, 
can be applied. 


4. Dynamic Stability Regions 


An analysis of the dynamic stability regions depending on the aerodynamic 
and mechanical parameters of shells consists in determining the boundary separa- 
ting the region in which s roots of (3.16) have negative real parts and from that 
where £,—s roots have positive real Darts, Sul, Da, e(Chw Poe [Ze 
To this end we insert w = iw in (3.16), where the real number @ varies within the 
limits — o<@ <o, and obtain by virtue of (3.16): 


(4.1) SEE) (io) = EEn0) +iEEn,0) = 0, 


where &, are arbitrary parameters of which the influence on the flutter range 
is the object of our investigation (€ = U, y =I, for instance). 
From the relation (4.1) we have the system of equations 


(4.2) E(é,©) =0, E(é,y,@) = 0. 
Eliminating w by means of the Fermat method, we obtain: 
(4.3) A(é,n) = 9 

and, for instance, 

(4.4) G(é,m) = 0. 


On the basis of (4.3), (4.4), and computing numerically &,7,0, we obtain gra- 
phically the boundary separating the regions where s roots of (3.16) have Rewa<0. 
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5. The Displacement Functions v,v,w, and the Flow Potential ® 


On the basis of the relations (2.1), (3.8), (3.2), the displacement of the middle 
surface u, v, w of the shell and the perturbation potential of the supersonic flow 
@® are expressed in the form 


f a—— DY CLM citarsneson, 
| y=] v! 
_ Ves) i(a,x-+np-+ ot) 
y = Gag eee NP+ Ot) | 
(5.1) =! 
footie 
ai 
<1 (wa U)o : ; 

0) — (o+4,U)o | mr el(a,xt+np+ ot) 

[ 2 v!m,I,(im,R) nit) 


The constants C, will be obtained from the normal system of equations (3.9) which 
has, as is known, ¢« independent solutions (¢ is bounded and positive). Similarly 
to the case (3.9) it can be shown that the non-homogeneous system of equations 
obtained from (3.9) is fully regular, therefore it has a unique set of solutions func- 


tions (assuming that C, is bounded). From the above it follows that all the series 
(5.1) are convergent. 
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DRGANIA POWLOKI O SKONCZONEJ DLUGOSCI 
PRZY NADDZWIEKOWYM PRZEPLYWIE WEWNETRZNYM 


W pracy autorzy rozwazaja problem drgan samowzbudzanych cienkiej powloki 
cylindrycznej o skonczonej diugosci przy naddzwiekowym przeplywie wewnetrznym. 
Zagadnienie jest rozwiazane przy warunkach brzegowych dla najezesciej spoty- 
kanych w praktyce odpowiednich urzadzen: jeden koniec zamocowany, a drugi 
swobodny (przy innych warunkach brzegowych zagadnienie mozna rozwiazac 
analogicznie). 

Rownania powtoki sa wyprowadzone wedltug teorii zagieciowe] z uwzglednieniem 
thumienia wewnetrznego zgodnie z modelem Voigta analogicznie jak w pracy [5]. 
Zachowane sa jednak pewne cztony w rownaniach r6wnowagi podobnie jak w [10]- 

Przeptyw wewnetrzny jest rozwazany w oparciu o zlinearyzowana teorie¢ nieusta- 
lonego przeptywu. Metody podane w tej pracy oparte na sprowadzeniu zagadnienia 
do nieskoniczonego uktadu réwnan algebraicznych pozwalaja na efektywne, nu- 
meryczne przeprowadzenie obliczen (miedzy innymi granic statecznosci dynamicz- 
nej) i na ewentualne zastosowanie maszyn matematycznych. 


Peswme 


KOJIEBAHHA OBOJIOUKM KOHEUYHONM JVJIMHBbI 
IIPH CBEPX3BYKOBOM BHYTPEHHEM TEUEHUU 


PaccmaTpuBaiorca camoBosbyxKqaomjue KoueOaHHA TOKO WMWIMHApHueckolt 
oO6ONOUKH KOHCUHO! TJIMHEI CO CBepX3BYKOBOM BHYTPeHHeM TeueHHen. Sayaua 
pemlaeTca Ip KpaeBbIxX yCNIOBMAX Id HaMOosIee uacro BCTpeualomyuxcaA, Ha 
IPaKTHKe, COOTBETCTBYIOWNMX YCTpOlicTBax, T.e€. Kora OAM KOHEI, SalljeMseH, 
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a pyro cBoOoqHEIa (pH Tpyrux KpacBbIx yCIOBHAX 3a]jauy MO%KHO PelIMTb 
aHaJIOruuHo). 

YpaBuenusd OOON0UKM BbIBOJATCA CormacHo Teopuu u3sru0a C yUcTOM BHYTpeH- 
Hero TemMMdbupoBanusa, cormacHo mMoyemm Dolirra, ananoruuHo padote [5]. Co- 
XPaHAIOTCA OHaKO HEKOTOPble UWICHbI B YpaBHeHUAX paBHOBeCHA MoOOOHO 
KaK B paOore [10]. 

BuyTpeHHee TeueHMe paccmaTpuBaeTCA, OCHOBbIBatch Ha JIMHeEapH30BaHHOw 
TeCOPHH HECTaljMOHapHOrO TeUcHHA. 

Meroybl, IpMBeCHHEIe B 9TOM paboTe, COCTOAT B cBeeHHM 3ayjaun K OecKO- 
He€YHOH cucTeme amreOpaMmuecKux ypaBHeHHii, MO3BOJIAIOT Ha e@eKTHBHO2, UHC- 
JICHHOS TIPOBE]CHHe PacueTOB (MEOKTY UPOUHM MpeseOB WuHAaMUYeCKOH ycTOli- 
UHBOCTH) MH BO3MO%KHOC MCHOUb30BaHHe MATCMATHUCCKUX CUCTHBIX MauIvH. 
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LINEARIZED SUPERSONIC FLOW PAST A VIBRATING SURFACE 
GF A BODY OF REVOLUTION 


ZBIGNIEW DZYGADLO (WARSAW) 
1. Introduction 


The aim of this paper is to solve the problem of linearized supersonic potential 
flow past a body of revolution, whose surface vibrates in a harmonic manner, and 
to determine the pressure of the gas acting on the vibrating surface. The solutions 
obtained will be applied in next papers to the investigation of self-excited 
vibrations of cylindrical and conical shells of finite length. 

The potential of a pulsating source in a homogeneous supersonic flow, which 
is known for very long and applied to many cases of flow, [1], [2], [3], can also 
be used as a basis for the construction of an axially symmetric flow potential past 
the vibrating surface of a body of revolution. 

Considering such sources distributed along the axis of the body W. H. DorRAN- 
cE, [4], and N. F. Krasnow, [5], obtained also the potential of distributed dipoles 
proportional to cos a (where a is the angular coordinate of the cylindrical reference 
frame in which the problem is considered), and applied it to the analysis of the 
flow past a vibrating rigid body of revolution. In the general case, the form of the 
potential of the unsteady supersonic flow past a pointed slender body was given 
in the first approximation by J. W. Mites, [3], valid only in a sufficiently small 
neighbourhood of the surface of the body. Due to the simplifications made, the 
only case considered with this approximation is that of perturbed transverse flow 
treated as incompressible, the axial variable and the time being treated as para- 
meters. 

The approximate form of the potential of a supersonic flow past a vibrating body 
of revolution of a considerable length was determined by G. I. Kopzon, [6]. The 
usefulness of the approximation of this seems to be somewhat doubtful if it is borne 
in mind that (as follows from the Eq. (9) of [6]) the rejected part tends to infinity 
as 1/y?, if y > 0, that is if the point on the surface of the body tends to the front 
edge and the remaining approximate expression for the potential is bounded in 
the neighbourhood of the point y = 0. 
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The Laplace transform of the potential of the supersonic flow past a vibrating 
cylinder has been given recently by M. Hoit and S. L. Strack, [7]. However, 
the effective determination of the original was not considered. 

In the present paper, the general form is determined of the potential of the ex- 
ternal supersonic flow past the vibrating surface of a body of revolution. Making 
use of the potential obtained, the solution of the problem of flow past a vibrating 
surface of revolution is reduced in the general case to the solution of a certain Vol- 
terra integral equation of the second kind. For a vibrating cylinder, the relation 
between the pressure and the normal component of the displacement is derived 
in the form of an asymptotic expansion in ascending powers of the inverse Mach 
number. From the expression thus obtained, we obtain, in the first approximation, 
an equation for the pressure in agreement with the so-called piston approximation. 
In contrast to the plane flow, where the error of the piston approximation is of the 
order 1/M? in relation to 1(M is the Mach number of the unperturbed flow), the 
error in the case of a cylinder will be of the order 1/M. In this connection, if self- 
excited vibrations of cylindrical shells are considered, further terms of the expan- 
sion should be used for greater accuracy, with moderate Mach numbers (MVM < 10). 

For pointed slender bodies the second linear approximation of the potential 
is given. It has an error of order double that of the first approximation, [3], takes 
into consideration the compressibility of the flow, and enables a more accurate 
description of the phase angle between the normal component of the displacement 
of points of the vibrating surface and the pressure, which play an essential role in 
the investigation of self-excited vibrations. The potential derived is used for the 
determination of the pressure on a vibrating cone. 


2. Fundamental Equations 
Let us consider an external supersonic flow past the vibrating surface of a body 
of revolution the OX-axis of which is parallel to the unperturbed flow of velocity 
Uy > a (Fig. 1). 


Bigsel 


The form of the body the surface of which is determined at rest, by the equation 
Yo = To(x), and the Mach number M of the unperturbed flow satisfy the conditions 
dro 


3 —.—— dr 
2.1 =e 2 0 
(2.1) | j <1 VM 1 <0 <ie 
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Assuming that the amplitudes of vibration of the points of the surface of the 
body are sufficiently small and bearing in mind (2.1) (if the number 6 is suff iciently 
less than 1), the problem may be solved on the grounds of the theory of linearized 
supersonic potential flow. 

Let us introduce dimensionless quantities. The coordinates x, r and the normal 
component of the displacement of the surface of the body W will be referred to L, 
the time ¢ to the quotient L/U, and the perturbing potential ® to the product 
U,L. For the reference lenght L, we shall choose the maximum radius of the body 
or its length depending on the necessity. , 

The dimensionless perturbing potential ®(x, r, a; t) satisfies the equation: 
oD PH 100 1 OG OD OD 


+2M? pee 


DD a i 
oe) Ox or? r or  r* 0a? Oxdt 


Paws 


where uw = V M1. 
The boundary condition of the Eq. (2.2) coupling the vibration of the surface 
with the flow has, bearing in mind (2.1), the form: 


aD dr, oD a. a 1 /dry\2 
oy [PE] (Sed ri] ren 


In view of the supersonic character of the flow, the potential © must also satisfy 
the condition 
(2.4) for x<pr, ®(x,r,a;t) = 0. 


The potential ® being thus obtained, the difference AP between the pressure Py in 
the unperturbed flow and the pressure P in the perturbed flow will be determined 


from the equation 
0b  o0@ 
Epp. Been Pate 
(2.5) NP AP, P ali ar 4 ml 
where Q, is the gas density in the unperturbed flow. 

In what follows we shall be concerned with the analysis of the vibration 
(2.6) W(x,a;t) = cos na ety, (x). 

In the connection the potential will be sought for in the form 
(2:7) P(x,r,a;t) = cos na exp[i(wt—fx) ]y,(x.2), 
where z= yr, B = M?w/M?—1, n=0,1,2,..., and the frequency of vibration 
@ may, in general, be complex. 

Knowledge of the solutions of the form (2.7) makes possible to study 
vibration of the form (2.6) and, solutions of more complicated cases may, by 
adding such expressions, be obtained. 

Substituting (2.6) and (2.7) in (2.2) to (2.4), and denoting 


M P | mead her RO Ie drone : 
Ca). 7 = > Wax) = 2 ot > | iol D (zs w,(X), 


3 Problemy drgan 
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the flow problem under consideration reduces to 


Op, Op 1 dy EM (he 
Sie et est Ms n ey = 0, 
Ce) Ox? O24 BeenOZ +(y T ze | 
Op, 1 dro Op, , iP Aro . 
2.9) | Z Se ae. 
OL) Vata xk en UeOs 5 
(2.10) for Xa.7  0,(x,2) = 0; 

The pressure difference determined by (2.5) may now be written in the form 
10) NP(X,a,231) = 0,U2 cos na emt A p, (x42), 
where the function Ap,(x,z) is expressed thus: 

-2.( OP o) 
=p —ibx n ; ; 
(2.12) Ap,(x,z) =e e iS | 


3. Solution of the Eq. (2.8) 


Let us proceed now to find the solution of (2.8) satisfying the condition (2.10). 
We apply the Laplace transformation to (2.8) in relation to x, and denote 


co 


Lo, (x,z) = i er 70, (e,zax == © (8.2). 
ry 


Assuming that ~,(x, z) is continuous for x >0 and z >0, and satisfies (2.10) and 
(2.8), we obtain: 


= = ; 2 
(3.1) Eh ets [se+9%4 | Pa= 0. 


The solution of the Eq. (3.1), the original function of which satisfies the condi- 
tions (2.10), has the form 


(3.2) Gn(8, 2) = A,(5)K,(z V2? +y?), 


where 4,(s) — is an arbitrary function, K, — a modified Bessel function of the 
second kind, the asymptotic appraisal of which for z—+ o is 


nom Eefisl 


The solution of (2.8), satisfying the condition (2.10), will be obtained by 
treating (3.2) as a convolution and writing its original function in the form: 


6.3) tal%52) = f AONAO—E, 2d, 
0 


where 


N,(x,z) = 27K, (z fs*@+9?). 
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The kernel N,,(x, z) of the expression (3.3) will be determined by applying the Efros 
theorem 


for x<z WN,(x,z)=0, 


(3.4) 


nl arch — 2) : 
LOL SZ oN (42) = = —yx | ng 2—B)x 


x2_72 


[ z 
ch(n arch x ch (” arch — 5 2) 


— ——— 
Vy 22 Vj e—8 ety | V @—z 


5K a, 


dé i Lov FP) Alyo—m] 


where J, and J, are Bessel functions of the first kind. For n = 0, the Eq. (3.4) may 
be transformed to give: 


0 forixy=z: 
(3.5) No(X,2) =) cos(yy/ x®—z?) 
Se 


Substituting (3.5) in (3.3), we obtain: 


fOmee=z. 


cos(yy (x—£)?—2?) dé. 
Va}? 


G.6) p(X. 2) = | Ad(é) 
Q 


Now, after substituting (3.6) in (2.7), we obtain the known potential of linearly 


distributed pulsating sources in a homogeneous supersonic flow [1]. 
In the case of n >0 we can obtain a solution of the Eq. (2.8) in the form (3.3) 


the kernel of which is of a much simpler form than that given by the Eq. (3.4). 
For this purpose, the following method will be used. 

There is a well known property of the Eq. (2.8), used in the theory of flow past 
bodies of revolution, consisting in that its solution for 1 = 1 can be obtained by 
differentiating with respect to z, the solution for n = 0. 

This property may be generalized to the case of any n = 1, 2, 3, ... For, 


eC? 1 0 (n—1)? Geo lee nn 
oo) I 6z? = Zz OZ z |= (aa 'z Oz a 
if 
(3.7.1) p= oe ss zon 


Applying n times the Eq. (3.7) we obtain 


q2 G2 1 @ 2 
(3.8) Lf i eel pee ae 


Oz te az OZ oz" ZO mee 


3* 
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where 
: Olcott tate. 3 1-2? 
— Sena acato ee ee d 
Cen) Pe Le ~éz z dz 2 oz al a 
From (3.8) it follows that if a sufficiently regular solution of (2.8) is known for 
n = 0, we can obtain the solution for n = 1, 2, 3, ... in the form: 
(3.9) PnlX» ) =o Lol; Z). 


Let us use (3.9) to determine a simpler form of the kernel of the potential ~,(x,z) 
than that given by (3.4). Let us assume (3.6) for the initial potential v(x, z). In 
view of the necessity of manifold differentiation we shall write it in the form 


x 
arch — = 


(3.10) (x, Z) at Ay(x—z chu) cos(yz shu) du 


obtained from (3.6) after changing the integration variable § = x—zchu. Sub- 
stituting (3.10), in (3.9) and assuming that A,(é) is differentiable a sufficient number 


of times and A,(0) = 4,(0) = ... = Ay” YO) = 0, we. obtain: 
ys arch = 
GlQe ox, 2) = L(x 2) = Al A (x—z ch uv) [B,, cos (yz sh u) ch nu+ 
k=0 
+C _sin(yz shu)shnu]|du 


where B,, and C,, are constants. 
Each term of the right-hand member of (3.11) identically satisfies (2.8). There- 
fore from (3.11), two different expressions for ,(x, z) are obtained: 


x 
arch — 
Z 


BA2)I" Oo (Kez) = J Gin(x—z ch u) cos (yz shu) chnu du = 


acall din(€) ln 5)" male (x—é,z)d— forn=0,1,2,... 


‘ V G-H-2 


arch ~ 
Z 


(3,13) (x, 2) ii Jon(x—z ch uw) sin (yz sh u)sh nu du = 


0 
I—2 


=f tn) nly oo ee Sas b)adt. for m= 102 oars 
( 


0 
where 


anloe Zz) = 5 lot P24 Py, 


fas, 2) = 5-10 +V/ 22x PAI, 
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4n(€) and qz,(é) are arbitrary functions. 
We can also construct the potential g(x, z), completing (3.13) for n= 0. 
It has the form 


ay 
arch— 
Zz 


(3213 1). 06322) = | Gao(x—z ch u) sin (yz shu)udu = 


sin [y/ (~—&?—22 Ep] eee ee : 
=f oo(5) (Gob “| 3 +7/ 5 g, 


The two forms (3.12) and (3.13) of the potential are equivalent. In what follows, 
one form or other will be used as required. 

For the applications, it will be necessary to know the Laplace transform of the 
potential (3.12). In view of the uniqueness of the problem the transform sought 
for may differ from (3.2) by a form of the function of the parameter s only. There- 
fore it can be determined by equating the transform of the asymptotic expression 
obtained from (3.12) for z—+0 with the corresponding asymptotic equation ob- 
tained from (3.2). In this way, we obtain for z—0 and n = 1, 2, 3, 


ae a(t) WV callie fanlx—&, 2)dE = 
0 


V @&—é — 


~ pe | dul eos ty OME, 


0 


gr-1 e 2(n=1) Re’ (s-tiy)" 
0p [ Mul e—O 008 yo ME = SON eee 
0 
and 7 
217-1)! A, ( 
A, ats) A. (zy Ase s?-y?) w ez eo : 
Hence 
i Nae Pa Re’ (s+iy)" 
A, ($8) = Iu )—73a53 (sh y2)2- 
and, finally, 
te ACE 
G1) 25 [ 1 Sa Fale b2t = 
: Re’ (s+iy)"” 


— Wan (s) (ep K,(zy s?+-y?) c 


where Re’ denotes the real part assuming that s is real. 
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4. Solution of the Flow Problem 


We shall discuss now a method for solving the problem of external flow. This 
consists in reducing the problem to the solution of a Volterra integral equation 
of the second kind. 

The solution may be obtained by using the potential w,(x, z) directly in the form 
(3.13). Making use of (3.12), we can also construct a potential with the kernel 
R,(x, z), bounded and non-zero for x = z and use it to solve the problem. 

In the first case, substituting (3.13) in the boundary condition (2.9), we obtain 


the equation 
X—Zo(x) 
M,, y=, Z9(x)] 


q2n dé = — ee td > 
TBP atay a Oe = — Kesey h) 


(4.1) 


where M(x, z) is a function bounded and non-zero for x = z; 


iC) nes x4 ; 20) ipy x? — —ain(= + -1)|l+ 


i age al In (24/2 -1) 
+> eS (y/28— 2) — sin. (yy x 2 | ! x as) Zz y 2 


V 2-2? fe dx 


sin (yy x22?) Zz aan dry 


M(x, 2) = fits af ae rere es Aa 3)| + 


see Re) Wesel | epmecen ies Che 2 a 
+) cos (yy'x zt )— (tae (: oF Pe > ene 2 


Substituting x, = x—z,(x) for x in (4.1), and applying the operator fe 


ym 


a9 


we obtain a Volterra equation of the first kind, which becomes, on differentiating, 
the sought-for equation of the second kind. 


The above method is convenient for bodies of simpler forms (cylinders, cones). 
For more composite forms another method is better suited. It consists in applying 
a potential with the kernel bounded for z >0 and non-zero for z = x. 

Such a potential, 

x—-Zz 
(4.2) Pals 2) = [ B(ER (XE, 2)dE 


10) 


{ 
: 
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can be constructed by means of (3.12). Indeed, assuming Gi.(S) = ly E, we obtain: 


in! is cos [yy/(x—£)?—z? ] 
(4.3) a a er 
0 


The function R,(x, z) is bounded for z >0 and non-zero for x = z, as can be seen 

by introducing in (4.3) a new integration variable by means of the relation 

€ = (x—z)o. 

cos {yy/(x—z)1—o) [x —o) +201 -+0)]} 
vod —o) yx —o)+z(1-+o) 


1 
(4.4) R(x, z)= call finlx(l1—o)+-z0,z]do, 
0 


4.5) RZ, Zz) = for x—>z. 


1 
do +e 
V 2z V o(1—o) V 27 
0 
Substituting (4.2) with the kernel (4.3) in the boundary condition (2.9), and 
bearing in mind (4.5), we obtain directly the Volterra equation of the second kind: 


(4.6) 8n(X1)— i (OH (X15 dE = hyo), 
0 
where 
rs V 2z1x)] : 
h, (4) an ee oft i dro{x(x,)] Wr [x(x,)], 
| pee ep de 
dx 
xf Ja OR, (x—E, z) 1 dro(x) OR, (x—E, 2) 
Ay(%, §) “a Lair 3(x) | z (OLD Ox 
a ianar) 
eax. 
me GND, R,(x | 
M eee 


and x(x,) is the inverse function of x, = x—zZ,(x). 
The solution of (4.6) will be obtained in the form: 


(4.7) Sn(X) = Iq) + ‘i H*(x,, Oh, (dé, 
i 0 


where 


Ay (X41, é) = Hm, aes Fas é) 
m=1 


and 


Tear) = : H,,(%1, DH am a(t, Oat. 
: 
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5. Determination of the Pressure Acting on a Vibrating Cylinder in a Supersonic Flow 


In the case of a cylinder we can, without solving equations of the form (4.1) 
or (4.6), obtain directly, in an explicit form, the relation between the pressure and 
the normal displacement component of points on the surface of the cylinder. 

Let us assume for the case now under consideration that the reference length 
is L = R, where R is the radius of the cylinder, and that the vibrating cylinder 
constitutes a portion of a body of revolution satisfying the conditions (2.1). The 
reference frame will be assumed according to Fig. 2. 


Fig. 2 


Let us apply to (2.12) the Laplace transformation in relation to x. In the relation 
thus obtained we substitute ¢,(s, z) in the form (3.2). Then, we obtain 


(5.1) Ap,(s, 2) = (s+iw)A,(s+iB)K,[zV (s+ip +7? J. 
Next, assuming that w,(x) is continuous for x >0 and w,(x)=0 for x <u 


(Fig. 2), and bearing in mind that now z,(x) = mw, we obtain, after transforming 
the boundary condition (2.9), the relation: 


(5.2) Tigh eee ee 
; py (s+iB+y?K, (uy (s-+iBP +7 ) 
Substituting (5.2)) in (5.1) on the surface of the cylinder (for z = jw), we obtain, 


(5.3) Bp, (s,) = — Se Sl 
pV (st+ipp+y Kiluy 6+ipy+y | 

The original of (5.3) can of course be sought-for by making use of the residues at 
the singular points and the branch points of the above relation. In this way an 
accurate equation is obtained. However, its form is useful for numerical compu- 
tation only and is inconvenient for other applications for instance for the analysis 
of self-excited vibration of cylindrical shells. 

A different form of the original, very convenient for applications, may be ob- 
tained by making use, for the inverse transformation of (5.3), of the asymptotic 
expansion of the functions K,(z) and K,(z) for large z. Such a procedure is justi- 
fied if it is assumed that the number M is sufficiently high. Then, the coefficient 


p= i. M?—1 appearing in the argument of K, and K/, may be treated as a large 
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parameter. It is true that this makes impossible the consideration of a certain range 
of M near to unity. That, however, is no serious limitation on applications of prac- 
tical importance. The asymptotic expansion K,(z) has, for large z the form [8]: 


4 m—1 k 
caf) Sapa 
k=0 


and we must have m >n—1/2 and 6| <1 if Rez>0, 
(4n?—1)(4n?—32).. .[4n®—(2k—1)?] 


@,0=1, @h= ae 


fork a1 2334... 


From (5.4), we obtain the following expansion of the function K,(z)/K/(z) if 
m >n-+1/2 


m—1 
(5.5) K,@) _ =| (2, kK), p(n, | 
k= 


K,@) Gee 2 
where (n,0)=1, @,1)=—1, 
k=1 
(1, k) = (0, k)—-(@, ky (n,j)@,k—jy for k=2,3,4,..,m—1, 


j=l 
and 


m—1 
(n, m) Te (n, m (n,m) 5 SNe ih) 
j=1 


and 


Gk ra (4n?-+-4k2—1)(n, k—1). 


On the basis of (5.5) we can write the asymptotic expansion of (5.3) for high 
Mach numbers: 


m—1 


a 2 (s+-iw)’w,(s) (n,k) | 
(5.6)  Ap,(s, “) uy (SLiBP Ly [ a [2uy/ (s+ipP+y?} 


: [Quy (s-+Hipe ry? 


Performing the inverse transformation of (5.6) term-by-term, we obtain(*) 


m—1 n, k) 
(5.7) Apy(.1)= orate +a) J i ee niches] S Ga seve 8+ 


(n,m) 


d. 
(Qu yn? Yn (X | ae 


() After the inverse transformation of (5.6), the reference frame is shifted to the right by 
4, so that the lower bound of the integral in (5.7) is zero, which is more convenient for the 
applications of the above relation. 


O(n, m) | 


+O 
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The function p,(x) have the form: 


! 2 
for k= wea) = 2 tyr t= oe (TY Kom), 


\ 


tos x ytL 
/ 1 2 
for k=29—1  @iyaG) = 2 14) 7s = i eS S,.40%), 


where y = 0,1, 2,3, ... and J, is the Bessel function of the first kind and the y-th 
order. . 

From the expansion obtained it follows that if the accurate expression for the 
function Ap,(x, ) is replaced by the partial sum of the first m terms of (5.7), the 
error will always be of an order higher by 1/u than the last of the retained terms. 
For m >n-+1/2 the absolute value of the error will be less than the first of the 
rejected terms. If MZ -—> oo the error tends to zero. For a fixed value of M there 
exists a number m, such that the absolute value of the error is minimum. 

Numerical computations have been performed for the quantity 
(n,m) 


6 V(X) 


ser 100% 
(n, k) 
a Yy(X) 

As (2u)k 


characterizing the convergence of the initial part of the expansion (5.7). The results 
are given in Table 1. The computations have been performed for n = 0; 2; 4, 


Table 1 M =2;4, x=0,5; 1;.1,5; 2; 3; 4, m=S and 

M\yw 1 0 > 4 @ = |. For w >1 the rate of convergence is 
05 0.9% higher and for w < 1 it decreases slightly. 

1 0.06% 11% - | 199) °_ From the computations and the form of 

os 1s He 79) ° 2 the expansion it follows that for M > 1,5 
2 389, o ri the value of the function Ap,(x, 4) can be 

1 : 2 ay determined with a high degree of accuracy 

15 56%, by retaining only a few terms (3 to 5) of 
—"“_ (5.7) provided that the following inequa- 


4 2 0.035% 0.56% Boys lities hold: 


3 0% Ee forn=0 x/M<1 
4 IAA 48% and for n = 1, Diy ee a WIM <li. 


If the ratio nx/M (for n #0) or x/M (for n = 0) increases above unity the rate 
of convergence of the initial part of the expansion decreases. For sufficiently large 


values of the ratio nx/M or x/M the expansion becomes divergent from the first 
term inclusive. 


The relation (5.7) can be represented in a still simpler form, very convenient 
for applications. For this purpose, let us observe that assuming for the value of 
the function Ap,(x,u) the partial sum of m terms of (5.7), we obtain Ap, with 
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an error of the order 0 (1/u™), which is equivalent to an error of the order 0 (1/M™), 
because = \/ M?—1, : 

Since the Mach number appears in terms of the retained partial sum, not only 
in the coefficient ~ but also in B and y, therefore, with no increase of the order 
of the error, we can obtain a simpler formula for Ap,(x,) if the m-th partial 
sum (5.7) is expanded in a series in 1/M, the retained terms being those with 
coefficients of lower order than 1/M”+!. Then, we obtain: 


°(z) 
Be = ee Fiow,o] ae 


z m—2 k 
1 Fae) — ame 1 . 1 
tae | |e Howl eee De bh x 2948+ (soma) 
0 j=0 


k=0 


‘ where the symbol c(m/2) denotes the integer part of m/2, and b, and b*‘;—the 
corresponding constants. For m < 5 they have the form: 
1 3 1 1 Wee cP 


bees), Diet b= 5 bn = > bi, = iw, bia = a OP URAL. 
St ee, to, 
b3, = — z jo" 1 @y? (25 | +E nt | 


Let us observe that using the first approximation of the expansion (5.8) 


(5.9) AP r(x, L) = = ee ions) 


we obtain a formula representing the pressure difference identical with the known 
equation used for the investigation of the self-excited vibration of plates and shells, 
which is the linear part of the so-called piston approximation to the pressure dif- 
ference in a supersonic flow. In contrast to the plane problem, where, using the 
relation (5.9), the error is of the order 1/M? in relation to unity, in the present 
case of a cylinder the error is of the order 1/M, because in the expansion (5.8) 
there appear all the powers of the inverse M number. Thus, from the form of 
the Eq. (5.8), it follows that if self-excited vibration of cylindrical shells is inves- 
tigated, further terms of the expansion should, for increased accuracy, be taken 


for moderate Mach numbers (M@<10). 
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6. Approximate Potential of Flow past the Vibrating Surface of a Pointed Slender Body. 
The Pressure on a Vibrating Cone 


We shall now derive an approximate expression for the potential of supersonic 
flow past the vibrating surface of a pointed slender body of revolution. This ap- 
proximation may be called the second linear approximation of the known theory 
of flow past pointed slender bodies, [3]. 


Let us assume now the reference length LZ to be equal to the length of the vibra- 
ting part of the body. Assuming that the condition (2.1) is satisfied for constant 
6 sufficiently less than unity we can determine the approximate potential @,(x, z) 
by retaining in the expansion of the accurate expression in a series in z, a few 
initial terms. Of course the equation thus obtained may be used only for the 


values 
Z=.2Z,(x)+ex, <where Wae<<l- 


To obtain the expansion required we shall use the transform ¢,(s, z) in the 
form (3.14). Retaining the first two terms of the series representing the function 
K,, we shall obtain: 


Sy eee ROC 00 Cat eee We 
GD ger) =e ae apenas +00] 


for, Wh 


~- 2g. eye 2 
(6.2) (s,z)= fale Era [I 7s Ly2) 4 o¢tina)| for n=2, 


= a (s) AY 2 a 
OD Se ge" Se f +24 y|2In (3 V4 7) ye i t 


+0(z* In ah for p41 


(6.4) Gls,2) = Gls) ae (#4 7) per Z (ot v[m(gvet)+e—]+ 


+0(24 In a} TOC 7) 0! 


where C = 0.5772 ... — is Euler’s constant. 


Performing the inverse transformation of (6.1) to (6.4) and introducing the 
notations 


x 


(6.5) Q,(x) = | ‘OOP *cosly(s—@)MiylE)d = 2-Mn— 1)! 2 AEING, ( 


for “1, 2) dase 


0 


Q(x) = qo(X), 
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we obtain the equations for the sought-for potential: 


6.6) g_(,2) = [- in Cag int) JOnc9-+0ce9| for n=3,4,5,.., 


ae 


61) g.2) =| a ("+ gue] [esc0-+0G"n 9 for n=2, 


72 


(G8) 94,2) = : if ) (in =— ; (r+ aes) ers (+ 


a\r 
+a] ih Q! (€)G(x— dé + 0(z4 In 2 for n=1, 
0 


Z a? 
COAG )=[In3 ! (5 ifr ira) [O00 


if 


2 qd? 7 
: [ | rites : i] [2sc@@¢ 48 4-02" n 2) for n=0, 
P.& 
0 


where 
Vx 


ae) = [Peas Inx. 


0 


For the derivation of (6.8) and (6.9), it was assumed that 
(6.10) Q;(0) = 2,0) = Q,(0) = 0. 


The potenial of the flow will be obtained in the first approximation, [3], by 
taking the first terms of (6.6) to (6.8) and the first term and also the first integral 
term of (6.9) and substituting them in (2.7). Although the first approximation is 
successfully used for the solution of problems of vibration of bodies as a whole, 
it is less suited for the analysis of local self-excited vibration of a covering (of 
the front portion of a body of revolution, for instance) in view of a relatively 
low order of the error and the incompressible character of the flow. The second 
approximation given above involves an error of order twice as high, takes the 
compressibility into consideration, and enables a more accurate description of 
the phase angle between the displacements of the points of the surface and the 
pressure, which plays an essential role in the investigation of self-excited vibra- 
tions. 

Let us apply the Eqs. (6.6) to (6.9) to the case of a vibrating cone of apex angle 
20. The potential ¢,(x, z) on its surface, that is for z = z9(x) = mu tg Ox, will be 
determined by means of the above equations with error of the order: 


O[(utg#)*], or Of(wtgd)* In (u tgd)]. 
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The function Ap,(x,), determined by (2.12) and characterizing the pressure 
difference in the flow takes now, on the conical surface, the form: 


(611) Aris 0) =e aq y* | (get) +7] [e+ jer 


n [\ 4 
foton=2,.3.4e 


(6.12) Ap,[x, Zo(x)] = : . + 5 ween? lm (+ tg o>] alles + i) + 


E Wied 2 d 7 ae n 
+) tie) oreo+ gat terex| (+18) +74] Ge) [ OF oO 


for 7. 


(6.13) Apolx, 21(x)] = =} nf 9 we tePx| +g ue tetoa" In[ > wtedx}—1] x 


Lléeofor emacs secrete 


x (at | [ es" @ore—oae for n=0, 
0 


where 
ca 


O10) = 


Q,(x), G*(x) = e-**G(x). 


From the boundary condition (2.9) we obtain, by rejecting the terms of orders 
higher than those appearing in the initial equations (6.6) to (6.9), the following 
equations. 

A. If, in addition to the conditions (2.1), we have: tg? > 1/u? 


~— 9) d 2 
14) {agente (at) +] bore = tgd eulraes Wal) 
for. n= 2, 3-4: ax. 
(6.15) se 2 tg? 9x] In z tg 0x) : d | ig) 4 || + 
: gh? tex In| 5H te Ox) +5 xt) HOF) + 


1 d ; ; 
poy te ol (7, + +9 if Or’ (6)G* (x—)dé = tg Ox? eas + ionc| 


+iow | 


format. 


(616) [1+ yet teton4] n(n tex] — [I (7. +18) +2" ] ose 
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1 7, 
+ Hite? onl (7, +) +7 ij 0% (OG*(x—2)dé = tex ee ane) +o] 
0 


forin==0. 
B. If the conditions (2.1) are the only conditions satisfied: 


(6.17) ie n— tet Ox 4 + aga | (i +16) +9"] lore = 


dw 
= te + — zie) wes nl) tow] LOle vimana See 


(6.18) fees are =p ted" nl -u tg os}: | as i) 7 llorey+ 


1 Tela : 
+ ute? o[(aet zy +2] i! O*'()G*(x—Bdé = 


tg o(1 +5 te? 0} Z Reg tion forn=1, 


1 d 1 1 1 d 
6.1 — to2 \ a erie 2 9, 4-2 oa ee a 
(6.19) : tg? 0x Ina tex} ae + a tg? 0x [in (qa tex} all ay + 


. : | * d i oa) . y : *, 7 
+18) +2*]l05() —|tet Ox 5 —pat tat ont (+8) 24] [ O50" 
0 


—&)dé = tg (1+ io) ke: o) ion] for n= 0. 


The conditions to be satisfied by the solutions Q*(x) of the Eqs. (6.14) to (6.19) 
follow from the consideration of the equation obtained by substituting the accurate 
expression of the potential (3.12) in the boundary condition (2.9) as applied to 
the cone. 

Thus it is seen that if 

w,(0) = 0 and w,(0) #0 


we have g,,(0) = 0 for every n and in order that g;,(0) = 0, it is necessary that 
w,(0) = wh(0) = 0. 
Bearing in mind (6.5) and (6.10) we can therefore write 
(6.20) OL\— 0, O)= 0770) =.= 07"0O) =0 forrrn=i1) 233s. 
assuming that w,(0) = 0 and w,(0)0 
and also 05(0) = Q3'0) = 9, 
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where we must have 
wo(0) = w5(0) = 0. 

In addition, in view of the uniqueness of the boundary condition with the accu- 
rate expression (3.12), we must have 
(6.21) Or(xy) =0, if ww= 

The condition (6.21) is in some cases necessary for the uniqueness of the solu- 
tions O*(x), because, in view of the supersonic character of the flow, we have only 
the conditions for x = 0, and this point is the singular point of the equations under 
consideration. 

We proceed now to solve the Eqs. (6.14) and (6.17). The casesn = 1 and n=0 
will be considered separately. 

For n = 2 we obtain from (6.14): 


1 | 
(6.22) OO’) = hes: ten oe sion]. 
From (6.17), after satisfying (6.20), we obtain: 
cn Dt PEEP dD 
(6.23) 010) = — Tee | rar tem [ae 
0 


where k = 2/tg? 3 
The Eq. (6.14) for n = 3, 4, 5,... may be solved by means of Bessel functions. 
The solution satisfying (6.20) and (6.21) has the form: 


2 dw 
(624) O*(x)= = Site = aVxe iBx) ¥, (yx) [ [=e Wn(S) nC) 


dw, (§) 


vas 2 ei, (yé)dé— J; wo | lar ion)! eB Y(vE)dé\, 


where J, and Y, are Bessel functions of the first and second kind, and 


acs 4 n(n—1) 1 
4 V aes Be 


For the values of A, for which the second integral in (6.24) becomes divergent, 
that integral should be replaced with the indeterminate one. 

The solution of the Eq. (6.17) for n = 3, 4, 5,... will be obtained by means 
of the confluent hypergeometric function ,F,(a, b, x). Satisfying (6.20) and (6.21), 
it takes the form: 


2(n—1)(2-+te29) 
(n—2)u?(b—1)tgd 


(6.25) On(x)= 


=D eerie 
Pan 2, F(a, b,—2iys) | & d 2 ei Btn X 
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X1F,@', 0-2 | e+ iom(@)| dt —x ies zat 2 Fa’, b',—2 iyx) x 


v8 


fee ay 1 DB +né nee (a, b, —2i | Ul 1@ Ww fal, 
0 


where 
ee (2-1) | 1 2(M—1)M 


~ (n—2)2 ae? ie (n—2)u? | 


SE eer aie lo in 2m ze 0) sterol or 12+ io] 


(n—2)u? ted 


4 1 ois 
b= tong) (n—1)(n—2yt(n-+-7 1840) +180] n+ ee: 


a’=1+a—b, b’ =2—b. 

For values of /+-b/2 for which the first integral in (6.25) becomes divergent, that 
integral should be replaced with an indeterminate one. The form (6.25) of solu- 
tion of the Eq. (6.17) is valid assuming that b40,+-1,+-2,.... For b =0,+1,+-2.,..., 
the solution will have a different form, which can easily be determined. 

On substituting (6.22) to (6.25) in (6.11) and then in (2.11), we shall obtain 
the sought-for pressure difference on the surface of the vibrating cone. 
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ZLINEARYZOWANY OPLYW NADDZWIEKOWY DRGAJACEJ POWIERZCHNI 
CIAL OBROTOWYCH 


W pracy wyznaczono postaé potencjatu naddzwiekowego oplywu zewnetrznego 
drgajacej harmonicznie powierzchni ciat obrotowych. W przypadku ogdlnym 
problem opltywu sprowadzono do réwnania catkowego Volterry drugiego rodzaju. 


4 Problemy drgan 
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Dla cylindra drgajacego w oplywie podano zaleznosé cisnienia od skladowej nor- 
malnej przemieszczenia w postaci asymptotycznego rozwiniecia wzgledem odwrot- 
nosci liczby Macha. Wyprowadzono drugie liniowe przyblizenie potencjalu opltywu 
ciat smuklych zaostrzonych i zastosowano dla wyznaczenia cisnienia na drgajacym 
stozku. 


Pesrme 


JIMHEAPUSOBAHHOE CBEPX3BYKOBOE OBTEKAHME IIOBEPXHOCTU TEIJI 
BPANIIEHWUA 


Onpenesaerca cbopma MoTeHIMasla CBEPX3BYKOBOPO Hapy2KHOrTO OOTEKAHHA KOJIe- 
Oarolllelica rapMOHHUCCKM MOBepXHOCTH TeJIa BpawjeHuA. B oOuyem cyuae sayaua 
CBOJHTCA K HHTerpaybHOMy ypaBHeHu1o Borstreppu Broporo poma. Ina yusmAypa 
KOJICOOMCrocd NPM OOTeKaAHHM aeTCAH 3aBUCHMOCTh TaBICHHA OT HOPMAJIBHO CO- 
CTABJIAIONIeH NepeMeMICHUA B BUJ,e ACHMMTOTHUECKOLO pasIOXKCHHA M0 OTHOLMICHHEO 
K oOparHoctu uncna Maxa. Brrpoyqutrca Bropoe JmmHeltHoe mpHOmmpKeHHe NOTeH- 
IMasla OOTEKAHHA TOHKHX 3a0CTPCHHBIX TeJI, KOTOPOe NPHMeHAeTCA [IA ompe- 
TeIeHvA WaBJIeHHA Ha KOJICOMOLIUMCAH KOHYCEe. 
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A TWO-TENSOR METHOD FOR INVESTIGATING NONLINEAR SYSTEMS 


WLADYSLAW BOGUSZ (CRACOW) 
1. Introduction 


Many technical problems are solved on the basis of the theory of ordinary 
differential equations. Such problems are, among others, those of linear elastic 
vibration, electrical vibration, stability of equilibrium of motion, orbital stability, 
engineering stability, and stability of transitory processes of automatic control. 
The solution of each of these problems requires a separate method. If a system of 
differential equations is linear with constant coefficients, the solution presents 
no essential difficulty. In the case of a system of nonlinear equations, or linear 
equations with variable coefficients, approximate or qualitative methods are most 
often used, based on topological properties of the solution of systems of ordinary 
differential equations. 

The behaviour of the solution of a system of two equations on the phase plane 
was studied by H. Porncarg, [1], by means of contactless curves. These are curves 
having no contact points with the integral curves determined by the system of 
equations. If a contactless curve is closed all its points are points of entry or points 
of issue. If such curves can be determined for a nonlinear system it is easy to 
describe the behaviour of the solution of this system. Also A. LAPUNOv, [2], uses 
surfaces of which the points are either entry or issue points. 

The problem of the behaviour of the solutions in the case where there are entry 
or issue points on closed surfaces was analysed by T. Wazewski in 1945, [3]. 
Using the notion of retract, WAZEWSKI arrived at a number of conclusions in 
the case where all the issue points on a closed surface are rigorous issue points. 

In practice, the determination of contactless curves or Lapunov surfaces whose 
points are only those of rigorous issue is, in most common cases, very complicated. 
In the present paper, a method is proposed for investigating the behaviour of the 
solutions of autonomic systems by means of spherical surfaces with their centres 
at the singular point. On the surface of such a sphere may appear entry, issue 
and slip points. The behaviour of the solution may be deduced from the pro- 
perties of two tensors, doubly covariant and contravariant of valence (0,2). The 
method is illustrated by a few examples. 


4* 
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2. Basic Assumptions and Dynamic Interpretation of the Method 


Let us consider an autonomic system 


(2.1) X = F(X), 


where X is a vector with coordinates (X,...,X,,) and F — a vector with coordinates 
(F,,...,F,). The assumptions are as follows. 

Assumption 1. The functions F; are of the C! class in an n-dimensional Euclidean 
space E,,. 

Assumption 2. The point (0,...,0) is an isolated singular point. 

Assumption 3. The Jacobian of the functions F; in the space E, is different from 


Zero: 
OF, 
aX, 


[8] 


Assumption 4. The transformation y = F(X) has its inverse in E,. The inverse 
transformation will be denoted by X¥ = ()y). 

Let us denote the partial derivatives 0F;/0X;, with the substitution for X, of the 
inverse transform X; = p(V1,---,Yn) by aj(V1,---, Yn), and consider the system of 
equations 


(2.2) y=aly)y, 


where a(y) is a nXn—matrix of elements a,,(J1,.-.,),). With these assumptions 
it will be shown in Sec. 3 that the study of the behaviour of the solutions of 
the system (2.1) reduces to that of the solution of the system (2.2). 

The form of the solutions of (2.2) in the phase space E,, depends on the pro- 
perties of two tensors: the doubly covariant tensor obtained from the scalar pro- 
duct of the distance vector @ with the coordinates (),,...,y,) and the velocity vector 
V with the coordinates (j,..., y,) and a doubly contravariant tensor constituting 
a bivector of the vectors V and p. 

The method for testing the solution of (2.2) will be interpreted for a system 
of three equations in the phase space £3. Let us consider the system of equations 


(2.3) XS ROGGE) peas ers) 


The system (2.3) may be interpreted as an equation of motion in the space E3 
of a material point p with unit mass. The trajectory of the point p passing through 
Po Will be denoted by f(po,t), and the velocity vector of p by V. 

The coordinates of the vector V are determined by the right-hand members 
of (2.3). The kinetic energy of p is: 


(2.4) B= OMAHA. 
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Let us differentiate both members of (2.3) with respect to time, substitute XY, = V3 
and apply the transformation and the notations (4). We obtain a system of equa- 
tions in the form (2.2), that is 

3 
(2.5) Vi= air J2» V3)V; (i= 2 3)5 


j=1 


the system (2.5) determines the motion of the point P with unit mass in the space 
E;. The trajectory of the point P passing through P, will be denoted by g(P,, 2), 
and the velocity vector by V (3,4, yo, y3) 

The distance vector of the point P from the origin 
has the coordinates (y,, ys, y3) and will be denoted 


Y3 


by e. ( 

Let the point p, be mapped (Assumption 4) into P). att 
The trajectory f(p,,t) is mapped then into g(P,,t), and - tn 
the kinetic energy (2.4) is half the square of the distance . } ; 
of P from the origin. The derivative of the kinetic Ay, : 
energy with respect to time is expressed by the scalar , 


product of @ and V 


d : j é 
(2.6) dt (E) = eV = ViVi tVoyetVsVs- 


Substituting for y, the right-hand members of the equations (2.5) we obtain a 
function of the variables y,, y2, ys denoted by @: 


3 
(2.7) DP(y1, Yo ys) pV = D PiIidys 


‘j=l 
1 : 
where the coefficients b;; = ay (a;;+-a;;) are functions of 1, ye, ys. 


If the function (2.7) has a constant sign along the trajectory g(Po, #4), the kinetic 
energy of p is a monotonic function along the trajectory f(po,t). In Sec. 3 it will 
be shown that in this case the point p tends to the singular point or its distance 
from the origin tends to infinity. 

Note. In the following the expression (a the trajectory tends to the singular 
point or to infinity) will mean that the distance of the point p on the trajectory 
from the origin tends to zero or to infinity. 

The sign of the function ® is decided upon by the matrix |fb,;|| which is the 
symmetrized Jacobian matrix of the functions F; with the inverse transform accord- 
ing to the Assumption 4. If the matrix ||/5,,|| is of a definite sign along the 
trajectory g(P),t) the function ® has a constant sign. It may happen that the 
matrix ||5,;|| is of a definite sign in the entire space E3. In this case ® has a con- 
stant sing along every trajectory g(P),t) and the form of these trajectories can 
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easily be determined. Every trajectory either approaches the singular point or is 
of increasing distance from it, thus tending to infinity. 
Such in this case is the geometrical interpretation of the form of the trajectory. 
Let us consider a sphere with the centre at the singular point and with any 
radius Q: 


(2.8) o? = y2+yit+y}. 


The velocity vector V, tangent to the trajectory g(Po,t), is determined at every 
point M of the spherical surface. 


h Ye 
Fig. 2 


The scalar product of the vector normal to the surface of the sphere at the point 
M and the vector V is the function ®. If at every point of the sphere (2.8) the 
function ® is negative, the point P moving on the trajectory g(P9, f) enters into 
the spherical region. Every point on the spherical surface (2.8) is an entrance 
point. If this property holds good for every @ that is if the function ® is negative 
in £3, every trajectory tends to the singular point. 

Similarly, if ® is positive, every point on the spherical surfaces (2.8) is an exit 
point, and the trajectories tend to infinity (Fig. 2). 

The case where the function ® may change sign in the space EF; is more in- 
volved. In this case, entrance and exit points may appear on a spherical surface 
and the kinetic energy (2.4) must be non-monotonic on certain trajectories. 

Here, three regions may be discerned in the space FE: the region w+ where the 
function @ is positive, the region w~ where it is negative and the boundary region 
«°. The points on a spherical surface (2.8) belonging to the region w° are points 
of internal or external slip. The point P) moving along the trajectory g(P,,t) may 
remain, starting from a certain time ¢, in the region w+ or w . Then, its distance 
from the origin increases or decreases. The point P may pass from wt to w~, and 
vice-versa, a finite or an infinite number of times. In these cases the form of the 
trajectory is tested by means of the moment of momentum of P about the origin. 
The vector of the moment of momentum is a vector product of the velocity vector 


OO 
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V and the distance vector e. The coordinates of the moment of momentum are 
expressed thus: 


(2.9) Pa SS ViVi Vi (1,of seords. 2508), 


The functions y;;, after replacing », with the right-hand members of (2.5), are 
functions of the coordinates (y1, yo, V3). 

Let us assume that all the functions y,; become zero in a certain region A. 
From the definition of the functions y,, it follows that in the region A the vector 
V has the direction of the vector e and an attraction or repulsion zone will ap- 
pear in the region A, depending on whether the vectors V and e have the same 
or opposite sense. Let us assume that the attraction zone is contained in the regions 
@ , where the function ® is negative and the repulsion zone — in the region w~, 
where @ is positive. If the regions A, w+ and w~ are detarmined, we can, as will 
be shown in the examples, determine the form of the trajectory g(P,, 1). It may 
happen that the trajectory g(P5, t) lies in the region A. Along such a trajectory 
the following equation is valid: 


(2.10) V =p, 


where the function 4 depends on the location of the point P on the trajectory. 
From the assumptions (2) and (4), it follows that the vectors V and p may become 
zero only at the singular point and, therefore, the function A(j, 2, y3) is con- 
tinuous and has a constant negative or positive sign. 

From (2.10) it follows that the trajectory g(Po, ¢) is, in this case, determined by 
the equations: 


t 
fron Yay ys)dt 
(2.11) Ve pes (i= 1, 2,3). 


The trajectory (2.11) may be entirely inside the region wt or w~. This depends 
on the sign of A(j1, yo, ys). 

Let us consider the other case, where the vectors of the moment of momentum 
(2.9) are different from zero in the space E3. The vectors V and p are in this 
case different from zero and determine the plane where the point P remains mo- 
ving along the trajectory g(P,, t). It will be shown in Sec. 3 that as between the 
functions e, ®, the following relation is valid: 


figs 
(2.12) OE aly €0) Be 


where gy is the modulus of the vector p at the initial time t = f, y — the modulus 
of the vector of the moment of momentum, g-the rotation angle of the vector e 
in the plane determined by V and p. The investigation of the form of the trajec- 
tory reduces in this case to that of the sign of the integral in the Eq. (2.12). The 
functions under the integration sign depend on g and ¢. In spite of this relation, 
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the estimate of the ratio of ® to w is often possible. In practice, it is more con- 
venient to examine the form of the trajectory g(P,, f), in the case under conside- 
ration, by projecting the point P moving along the trajectory g(Po, t) on the 
planes (),, ye), (3), (Vs, Wi), the form of the trajectory of P being deduced 
from that of the projections. 


3. Mathematical Justification of the Method 


Let the autonomic system be described by the system of ordinary differential 


equations 
(3.1) X = F(X), 


where X= (X,,...,X,), & = (i;...36,)) a5 in sec. 1. Let-us make’ the. Assumptions 
1, 2, 3, 4 of the preceding section and consider the system of equations 


(3.2) y =aly)y, 


where a(y) is an nXn—matrix. Let us denote, as before, the trajectory passing 
through the point p,cE, at the time t = ft) by f(pp,t), and that passing through 
P cE, by g(Po,t). The point py is mapped, by the transformation defined in the 
Assumption 4, into Py. On the basis of this assumption, we can write: 
(3.3) &(Po t) =F [fo] and  f(Po, t) = lg(Po. t)].- 
We proceed now to prove the following theorems: 

Theorem I. If the trajectory f(po, t) is bounded, the trajectory g(P,, ft) is also 
bounded. 

Theorem II. If the trajectory f(p,, t) tends to the singular point, the trajectory 
g(P,,t) also tends to the singular point. 

Theorem III. If the trajectory f(pp, t) tends to infinity, the trajectory g(P,, 1) 
also tends to infinity. 

Proof of Theorem I. Let us assume that the trajectory f(po, t) is bounded —that 
s the distance of p on f(po, t), from the origin is bounded. 

It follows that there exists an r>0 satisfying the inequality 


(3.4) 0[0.f(Po: )]<r for te [to, 00 |, 
where e[0,f(po, 1] = De and (Xj,...,x,) are the coordinate of p. 
i=l 


By virtue of the Assumption 1, the functions F;(X,,...,X,) are continuous in 
E, and bounded in the closed region (3.4) Therefore, there exists an R>O such 
that 


(3.5) Fif(Po. DICR, 


if f(po, t) lies in the region (3.4). From (3.5) and (3.3) it follows that g(Py, t) is 
bounded. 
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Proof of Theorem II. Let us assume that the trajectory f(p,, t) tends to the 
singular point—that is the distance of {the point p on the trajectory f(py, 0) 
tends to zero. From the definition of the distance of the point p on the trajectory 
it follows that all the coordinates X,; tend to zero. From the continuity of the 
function F;,(Xj,...,.X,) and from the Assumption 4, it follows that F(X) tends to 
F(O) and F(0) = 0. By virtue of (3.3), the trajectory g(P,, f) tends to the singular 
point. The inverse theorems of I and II are also valid and the proofs are similar. 

Indirect Proof of the Theorem III. Let us assume that the trajectory f(pp, f) 
tends to infinity and that the trajectory g(Po, f) is bounded. Bearing in mind the 
inverse theorem of I, we reach a contradiction, because, by virtue of (3.3), the 
trajectory g(P», t) becomes the trajectory f(pp, f). 

With the theorems just proved, the test of the solutions of the system of equa- 
tions (3.1) for boundedness and for the quality of tending to zero, reduces to the 
investigation of the form of the solutions of the system of equations (3.2). 

3.1. The doubly Covariant tensor and the sufficient condition for the straightness of the 
system (3.2). Let us construct a symmetric Jacobian matrix by adding to the 
matrix a(y) the transposed matrix a7(y) and dividing by 2: 


(3.6) bo) 5 tal) +a70)) = |5,||. 


The following additional assumption will be made for the matrix b(y) 
(3.7) || 5;;|| 4 0 in the space E,. 


Let us consider a point P on the trajectory g(Pp, f). The vector of the distance 
e of the point P from the singular point has the coordinates ()j,..., ,), and the 
vector of the velocity V has the coordinates (j,,..., y,). The scalar product of e 
and V is 


@ — N'y.j, = eV. 
(3.8) 2 e 


1=1 

Replacing j, with the right-hand members of (3.2), the function (3.8) is expres- 
sed by 
(3.9) D(V15-05 Va) = bud; 

i,j=1 

Let us consider an arbitrary point P, in the space E,, and calculate the coeffi- 
cients b,, at the point Po: 
(3.10) Bi Oyses Yr) sae bi. 


Substituting bf, in (2.12), we obtain a homogeneous quadratic form 


i. 


(S11) D(),...; Vn) 7a, DS Ii 


i, j=l 
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It is known that the coefficients of a homogeneous quadratic form are trans- 
formed as the coordinates of a doubly covariant tensor. It follows that the geometric 
entity (3.11) is a doubly convariant tensor. At every point of the space E,, we 
can construct a geometric entity (3.11), thus determining a tensor field in E£,. 
Let us consider the case where the matrix b(y) is determined at every point of 
the space E,. Let us suppose that the matrix b(y) is positively definite at every point 
of the space E,. The right-hand member of (3.9) can, in this case, be reduced to 
a sum of squares of which the coefficients are positive. The function D(),...., y,) 
has the plus sign in the entire space E, except the singular point where it is zero. 
It will be shown that in this case the system (3.1) is straightened. To this end 
the Barbaskin theorem, [16], will be used. If in an open region contained in E, 
there exists a function u(x,,..., X,) having constant partial derivatives, and if there 
exists a number k?>0 such that the inequality 


n 
Ou 
Ox; 
i=1 


(3.12) F>R, 


holds, the system (3.1) is straightened. In our case, the function u(x,..., x,) will 
be determined by the equation: 


3 i 
(3.13) U(X4,...5 Xn) = — [FG Xn)? C 
i=1 
It can easily be verified that 
(3.14) 7 Ou ee 
. Ox; Po Sa Vi gerry ay 


i=1 


on substituting yi = F,, Let us take any constants 0,>0,>0, and determine the 
region G, by means of the inequality 


(3.15) 0<e< D> <a: 


i==1 
In the region to, the function ® is continuous, positive and has a positive mi- 
nimum which will be denoted by k?>0: 
(3.16) P(),...5 YJ2k* >0. 
The inequality (3.13) holds also in the open region G,. The assumption 4 


makes the region G, into the region G, open and contained in E,. On the basis 
of (3.16) and (3.14) the following inequality holds in the open region G,: 


(3.17) s Loa tS) 
Ox; 
it 
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Hence, on the basis of the theorem just cited, the system (3.1) is straightened. It 
will be shown, in’addition, that the trajectories g(Po, t) tend to infinity. The square 
of the distance of any point P on the trajectory from the origin is: 


(3.18) o* = S yj. 
It can easily be seen that: 
d 2 
(3.19) <= WD pyyeory Ve 


The function O(),,..., y,) is positive in E,, therefore the function 9? is an increasing 
function. Let us fix t = t, and determine the value of the function 02(¢) at t, by 
o;. Let us take any ¢>¢) and apply the theorem on the mean value of the function 
e7(t) of one variable: 


0? (t)—e7 do* 
5; = 2 2 
(3.20) are Fr o?(t) >ei>0, 


where t*e[t,, ¢]. 

If ¢ tends to infinity, the denominator of the fraction (3.20) tends to infinity. 
It will be shown that the numerator tends also to infinity—that is 0?(t) increases 
indefinitely. Suppose that o?(t) is bounded. The fraction (3.20) tends to zero, 
and the derivative of 0?(t) tends also to zero. On the basis of (3.19), the function 
®[y,,(t),..-, ¥,(t)] would tend to zero, but the function ® is continuous in £, and 
becomes zero only at the singular point. Hence it is concluded that y,(¢) tends 
to zero and, by virtue of (3.18), 0?(t) tends also to zero, which is impossible in view 
of (3.20). The other case, where the matrix b(y) is negatively defined, may be located 
in an analogous manner. 

In this case, the system (3.1) is’also straightened, and the trajectories tend to the 
singular point. From the above reasoning, we have the 

Theorem IV. If the matrix b(y) of the covariant coordinates of the tensor (3.11) 
has, at every point of the space E, a definite sign, the system (3.1) is straightened 
and the trajectories of the system tend to infinity if the sign of the matrix is po- 
sitive, and to the singular point if is negative. 

The geometrical interpretation of the doubly covariant tensor (3.11) is such. 
Let us consider an arbitrary point P, in E,, and calculate at this point the coeffi- 
cients b,, as before. Let us substitute the coordinates of P) in the right-hand mem- 
ber of (3.11). We obtain a geometrical entity which is an ellipsoid described by 
the equation 


(3.21) DOV?, oon Yn) = YS" Bi Vid; 

i, j=1 
For linear systems, the coefficients b,;; are constants independent of the point Po. 
The ellipsoid (3.21) has constant principal directions. It is concluded that linear 
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vibration is such that only the length of the axes of the ellipsoid (3.21) varies, the 
principal directions remaining unchanged. For nonlinear vibrations, the principal 
directions of the ellipsoid vary with the point Py. 

Let us consider now the case where the matrix b(y) is not of definite sign at some 
points of the space E, or in the entire space E,. The form of the solutions of the 
system of equations (3.2) is tested by means of a doubly contravariant tensor. 


3.2. The doubly contravariant tensor. Let us form of the vectors e and V, the bivector 
w. The coordinates of the bivector are determined by the equation 


O22) Viz = ViVi: (ay =, ey A, i egy: 
or can be expressed as the minors of the matrix 
(3.23) | War ecee ys | 

pees 


The bivector of which the coordinates are determined by (3.22), is called a doubly 
contravariant tensor. The condition of the bivector becoming zero, is necessary 
and sufficient for the vectors @ and V to be parallel. In this case the terms in the 
first row of the matrix (3.2) are proportionate to those in the second row. Let us 
assume that the bivector w becomes zero on the trajectory g(Po, t). The vectors 
e and V are parallel on this trajectory, and from the parallel condition we obtain 
the relation between the coordinates of these vectors 


(3.24) Vi = Mi, 


where A is a function of the location of P on the trajectory g(Pp, 7). 
From Eq. (3.24) we obtain: 


t 


i Adt 


(3.25) yi = yee” 


Let us examine the conditions of existence of a function A satisfying (3.24). Sub- 
stituting the right-hand members of (3.2) in (3.24) for »,, we obtain the system of 
equations: 

(Qi — Ait QyeVet + Qin Vn= 0, 
(372.6 Ne iii Dig MUR asd aig cia Soothe ens Eo) 
An Vit AyeVet +++ Gin —AVn = 9, 


where a,; are functions of the variables (j,, ..., y,). The condition necessary for 
the existence of the function / is that the determinant 

ay—A eee Ayn 
B2))..= aoe." © Oise =O; 


Ant ++ Any — A ! 


becomes zero. From the Eq. (3.27), it follows that the function A (if it exists) must 
be different from zero. The free term in the Eq. (3.27) is the Jacobian (c), of which 
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it has been assumed that it is different from zero in E,. It follows that if there exists 
a function A, satisfying (3.27) and (3.26), the trajectory (3.25) tends to the sin- 
gular point for negative A and to infinity for positive 2. 

It can easily be shown that the function @ has, along the trajectory (3.25). a con- 
stant sign identical with that of 4. For this it suffices to substitute (3.24) in (3.8). 

Let us consider now the case where the bivector w is different from zero in the 
space E,, except at the singular point where it is zero. In this case, the bivector 
determines, at every point P of the trajectory, a two-dimensional hyperplane in 
the space E,, passing through the singular point. In this plane lie the vectors V and 
erie 3). 


higas 


The vector V can be resolved in this plane into two vectors: V, normal to ep and 
V,, whose direction coincides with that of p. 

Let us introduce on this plane an orthogonal reference frame, determined by 
the vectors e and V,, and having its origin at. The coordinate of the vector V, is 
always positive. The coordinate of the vector V, is positive or negative depending 
on whether the vectors V, and e are of the same or opposite sense, and may be 
determined by the equation. 


(3.28) Vp= = 


The elementary displacement of the end of the vector p in the direction of V may 
be represented as a geometric sum of the elementary displacements ds, and ds, 
along the axes of the reference frame assumed. The coordinate of the vector ds, 
will be computed by means of the equation 


(3.29) ds, =odp, or ds, =V,dt, 
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where dp is the elementary rotation angle of the vector e about the singular point. 
From (3.29) we have: 


(3.30) odp =V,dt, or ii z 


“dt 0 
It can easily be verified that V,0 = |p| and, on substituting in (2.35), we obtain: 


d 
(3.31) = a 


From (3.31) it follows that the angle y is a monotonic function if the bivector 
is different from zero outside the singular point. Let us compute the ratio of the 
coordinates V, and V,: 


3.32 _& — 
Cee V,  — dodt 

Bearing in mind the relations 

_ do _ 2 4p 
we obtain from (3.32) 
do/dt @ 

3.34 = : 
eee aldpidty — [| 

Rearranging (2.39), we obtain the differential equation 
(3.35) 1 do @® dp 


g dt |p| dt’ 
of which the solution is 


7) 
@ 
lor” 


(3.36) © = @e° : 

where @ is the modulus of the initial vector gy and the angle m is measured from 
Qo. The functions @ and | q| in the Eq. (3.36) depend on the coordinates of the point 
P on the trajectory g(Po, t). In many cases, an appraisal of the ratio of these func- 
tions is possible, and (3.36) may be used for testing the boundedness or the perio- 
dicity of the solutions of the system of equations (3.2). 

The following conclusions may be drawn from the above reasoning, concer- 
ning the form of the solutions of the system (3.2). In the case where the matrix 
b(y) is indeterminate at some points or in the entire space E,, we have: 

Conclusion I. If the bivector w becomes zero over a certain region, this region 
constitutes an attraction or repulsion zone depending on whether it is contained 
in the region ® < 0 or ® +0. The trajectory contained in the region |p| = 0 has 
the form (2.30). 

Conclusion IT. If the bivector p becomes zero at the singular point, only, the 
investigation of the form of the solutions reduces to that of the integral in (3.36). 
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In practice, it is more convenient, in that case, to examine the form of the tra- 
jectory of the point P by projecting P on the coordinate planes. From the form 
of the trajectories of the projections, some conclusions may be drawn on the spatial 
form of the trajectory. 


4. Applications 


Example 1. Vibration of a dynamic damper. The equation of motion of the 
dynamic damper represented at Fig. 4. with nonlinear damping characteristic can 
be written in the form 
(4.1) x+x+ax-+bx?+ cx = 0, 
where a, b, c are constant parameters of the ysstem 
and a >0. 

Substituting x = x, x =X, we obtain the system 
of equations: 


(4.2) a = Xa Xs —= —X,—ax,—bx3—cx3. 
Denoting yy =X, yo = —X,—ax,—bxi—cx3, we 
obtain the system of equations: Fig. 4 


(4.3) i=Ye Yo = —Yi—(at+-2by, + 3cyp)yo. 
The Jacobian of the system (4.2) is equal to unity. Let us calculate the function 
(4.4) D = —yx(a+ 2by,+3cy)). 


If b? < 3ac, the function ® is non-positive in the entire plane (j4, 2). By virtue 
of the Theorem I, every trajectory of the system (4.3) tends to the singular point. 
Let us investigate the coordinates of the tensor . Two coordinates of the tensor 
are zero, the third is expressed by: 


(4.5) y = —y2—yyyo(at2by,+ 3cy})—y3- 
If 
(4.6) a+2by,+3cyz < 4, 


the function yw is negatiwe. 
The inequality (4.6) holds, if 5? > 3c(a—2) and y, satisfies the inequality: 


b+y/ b?—3c(a—2 b—y b?—3c(a—2) 
asl a ) nc ——- 


(4.7) 


Outside the interval (4.7), the function y becomes zero on the curves: 


(48) 9 nt 22 @42by +30) FF V @t2by. + 3e7—4 = 0. 


A field of linear elements is distributed along the curves (4.8) thus constituting 
an attraction zone. 
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The form of the solutions is illustrated in Fig. 5. 

Let us consider now the case of b2—> 3ac. In this case the form of the solutions 
for c >0 and c < 0 is different. For c > 0 the function ® is positive, if y, satisfies 
the inequality: 


b+y b—3ac b— B® —3ac 


(4.9) Soy = 


BG 


Yoh 


OG 
¥-0 
Yo 
| 
BA 
ay $>0| | $<0 
aN | 
— 
¥=0 | 
TTT a 4 
lp eee 
—— 
Fig. 6 


The function y becomes zero on the curves (4.8), an attraction zone being formed 
as before. The form of the solutions is shown at Fig. 6. 

For c <0 the function © is positive outside the interval (4.9). The function 
y becomes zero on the curves (4.8) which lie in the regionD > 0, A repulsion zone 
is formed. The form of the solutions is shown at Fig. 7. 
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The integral curves of the system (4.3) in Fig. 5, 6, 7 are drawn in the following 
manner. We choose an arbitrary point P and draw the vector p. The rotation of 
e is counterclockwise in regions where y is positive. In regions where yp is nega- 
tive, e rotates in the clockwise sense. Drawing the integral curves in regions where 
® is negative, the distance of P from the origin is reduced and in regions, where 
® is positive it is increased. The form of the integral curves in the (—) vicinity 
of repulsion zones and attraction is obvious. 


Fig. 7 


From the integral curves of Figs. 5, 6, 7, it follows that the equilibrium is asymp- 
totically stable. The asymptotic stability in Fig. 5 and 6 takes place in the entire 
plane and in Fig. 7 at some places only. The form of the solutions of the system 
(4.2) is analogous, on the basis of the theorems 1, 2 and 3 above. The present 
example has been solved by the method described. However, the passage from 
the system (4.2) to (4.3) can be avoided, the system (4.2) being written thus: 


(4.10) x1 = X95 Ne = Xy X,(a | bx. Cx 


This system is of the type (4.3), and can be treated directly by the method de- 
scribed. The passage from the system of equations (3.1) to (3.2) is necessary only 
when (3.1) cannot be written in the form (3.2). 

Example 2. Let us apply our method to a self-excited system governed by the 
Van der Pol equation: 


(4.11) y—el—y)yty=0, e€ >0. 
Substituting y = y,, y = 2, we obtain the system of equations 
(4.12) Vi=Yo  Yo= el —Yi)V2—-J1- 


The system (4.12) is of the type (3.2) and there is no need to use the transformation 
from the assumption 4. Let us calculate the function @: 


(4.13) @® = cyz(1—y?). 


5 Problemy drgan 
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The function ® is non-negative in the zone —1 <ly, < 1 and y, is arbitrary. Since 
in this region lies a singular point, the equilibrium is unstable. Let us investigate 
the non-zero coordinate of the tensor wp: 


(4.14) y= —yi+yy2e(1—y)— J. 


The function y becomes zero on the curves described by the equations 


E —— 
(4.15) n—- SI —W+ GV ed —y-4 = 0. 


The curves (4.15) lie in the zones determined by the inequalities 


2 2 
(4.16) yay 4d. or n<-Wf 142: any ye 
(4.17) = oe ae joss any 2. 


There is an attraction zone along the curves (4.15), lying in the zone (4.16) 
and a repulsion zone along those in the zone (4.17). The form of the solutions is 
shown in Fig. 8. 


Fig. 8 


From Fig. 8 it follows that the trajectory starting from the point P, bounds 
a region into which every trajectory enters. On the other hand, the circle K bounds 
the region from which every trajectory exits. In this way the region of the inter- 
mediate process is determined. In the case where « is less than 2 the zone (4.17) 
does not exist and the curves (4.15) are the only remaining in the zone (4.16). 
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Example 3. A centrifugal governor of direct action, with non-linear damping, 
Fig. 9. 


The equation of motion of the governor and the engine can be written thus. 
(4.18) To=—p 


a 


— equation of motion of the engine, 
Trut+T,t+ke+he+ou =o — equation of motion of the governor. 


T, — acceleration period of the engine, 7, — natural vibration period of the 
governor multiplied by 2x, 7, — time constant of the damper, k, h — the 
parameters of the nonlinear characteristic of the viscosity damper, 6 — the 
sensitivity factor or the governor, y — variation of angular velocity, 4 — the 
displacement of the governor sheath. 


Let us introduce the notations 


k h 3 fT? Met, 
ery is P=. sat) os, B= a oe 


and the transformation 


With these notations, and applying the transformation the system (4.18) takes 
the form 


(4.19) 9=—x, x = g—Ax—Bx—aB*x*— BB x’. 

Substituting x = x1, X = Xy, P==X3 we obtain a system of three equations of 
the first order: 
oa = — (B+0B?x,+ BB°xi)x,;—Ax2+X3, 


Xe = X1> X3 — — Xo. 


(4.20) 


The system (4.20) has been written in the form (3.2). The investigation of the 
behaviour of the solutions will be undertaken assuming: 


(4.21) B?>4A, @—48<0, fB>0. 


5* 
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Let us denote the quadratic trinomial in the first equation by 
(4.22) W = B+aB?x,+fB°x%. 
With the assumption (4.21), the trinomial W is always positive. Let us investigate 
the motion of the projection of the point P moving along the trajectory, on the 
(x1, X2) plane, (Fig. 10). 


X3 
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SY ee 
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xy 12 


Fig. 10 


The scalar product of the vectors of velocity of the projection of V,, and the 
distance p,., is determined by the function ®,,, and the vector product of these 
vectors by the coordinate of the tensor w: 

(4.23) hee = —Wxi—x,x,(A—1) + Xe; 
Wie = XP + WyxXg+ Ax—X9X3. 


For every x, the function ®,, is zero on the lines 
(4.24) Xy = 0, Wx,+x,(A—1)—x, — 0. 


The regions in Fig. 11 where the function ®,, is positive are shaded vertically. Let 
us reduce the function y,, to the sum of the squares 


bea Wx,\? , 44A—W? Drew x 
(4.25) Vig (x 5) + 4 (= rer == 4A—Ww2- 


From the assumption (4.21), it follows that W* > 4A. Setting y,. equal to zero 
we obtain a non-parameter family of curves in the (x, x.)-plane. The variable 
x3 1s considered to be a parameter. The behaviour of these curves for a positive 
fixed x3 is represented in Fig. 11 a and for a negative x, in Fig. 115. Attraction 
and repulsion zones are formed along the curves y,, = 0. The regions where 74, 
is negative are shaded horizontally. The projection of every trajectory on a plane 


parallel to the (x,, x.)-plane, tends to the point Q with the coordinates (0 “a »} : 


At the point Q, the trajectories intersect at right angles the planes parallel to the 
(x1, X2)-plane. Let us investigate the motion of the point Q. To this end let us pro- 
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ject the point P lying on the trajectory on the plane normal to the (x, x,)-plane 
— the (x3, x,)-plane for instance. Let us calculate the function ®3,: 


(4.26) D,, = — Wxi— Ax 1X.—X9X3. 


a) Ye | 


Fig. 11 


Y4 


Fig. 12 


If the point P lying on the trajectory tends to QO, x, tends to zero, x, tends to 
x,/A, and the function ®,, tends to —x2/A. This means that the distance of Q from 
the origin decreases and x, tends to zero. The form of the solutions of the system 
(4.20) is represented at Fig. 12. The equilibrium is asymptotically stable in the 


entire space. 
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Streszczenie 


DWUTENSO ROWA METODA BADANIA UKLADOW NIELINIOWYCH 


W pracy opisana jest metoda badania przebiegu rozwiazan ukladéw nieliniowych 
autonomicznych w przestrzeni fazowej. Wykazano, ze przebieg rozwiazan zalezy 
od wlasnosci dwoch tensoréw, z ktérych jeden powstaje z iloczynu skalarnego wek- 
torow predkosci i odlegtosci punktu w przestrzeni fazowej, a drugi jest biwektorem 
prostym tych wektoréw. 

Metoda moze byé zastosowana do wyznaczania obszaréw statecznosci i do ba- 


dania charakteru punktu osobliwego. W pracy objasniono metode na trzech przy- 
ktadach. 
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Peswme 


UBYXTEH30PHbIM METOJL UCCJIEMOBAHUA HEJIMHEMHYX CUCTEM 


JlaetcA mMeTOA uccieqoBaHuA Mmporecca pewleHuii HeIMHeMHbIX ABTOHOMHBIX 
CHCTeM B (ha30Bom ImpoctpaHcTBe. JloKa3brBaeTcA, UTO Mporecc pelleHuit 3aBHCHT 
OT CBOHCTB J[ByX TCH30POB, H3 KOTOPbIX Me€PBbIM ABJIACTCH MpOUSBeeHHeM CKa- 
JIAPHOFO BeEKTOPA CKOPOCTH HM PaCCTOAHUA TOUKH B (ba30BOM MpOcTpaHCTBe, a BTO- 
pol MmpocTbim OuUBEKTOPOM STHX BeEKTOPOB. 

OTOT MCTO, MOET ObITh IPHMeHeH IA OMpenesenua oOnacTeH ycrohunBocTu 
M [IA UcceqOoBaHHA XapakTepa ocobol TouKu. IIpuBeneHHEIi meToy, wWWOCTpH- 


pyeTca TpemaA nmpHmepamu. 
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NONLINEAR THEORY OF A GENERATOR WITH A LONG LINE EXCITED 
BY A NEGATIVE RESISTANCE 


STEFAN HAHN (WARSAW) 


1. Introduction 


It is known that the amplitude of vibration in generating systems depends on 
their nonlinear properties. At the same time, due to the non-linearity of the system, 
the frequency of the generator depends on the amplitude. In practice, nonlinear 
influences on the vibration frequency can often be disregarded in relation to the 
linear factors (thermal, structural etc). However, even in this case the accurate 
knowledge of the nonlinear effects is essential, for instance, for the problem of 
amplitude stability. As regards generators with a long line, in spite of the fact that 
they are widely used in engineering — micro-wave engineering in particular 
the literature devoted to appropriate nonlinear phenomena is very scarce. Of im- 
portance is the paper by GroszKowskI, published in 1938, [1], in which the method 
of harmonics proposed by that author is applied to a generator constituted by 
a negative resistance exciting a long line. In particular, it was shown by GRosz- 
KOWSKI that in a generator with a long line the influence of nonlinear effects on 
the generator frequency may be considerably smaller than in generators with cir- 
cuits with concentrated constants. The work by W. Majewsk1, [4], should also 
be mentioned, in which a long line generator is considered, the method of diffe- 
rential equations being used(’). 

The nonlinear theory of a longline generator presented in the present paper is 
based on the method of the third harmonic proposed by J. GroszkowskI, [2], and 
used also in the former paper of the author, [3]. This theory should be considered, 
in principle, to be qualitative, explaining the mechanism of the physical phenomena, 
although quantitative results can also be obtained in certain particular cases. 


(*) This paper was published in 1938 in limited number of copies and is little known to the 
younger generation of scientific workers. 
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2. Admittance of the Negative Resistance 


Let us consider the generating system shown at Fig. 1. This system is compo- 
sed of a long line, short circuited at the end, to the input terminals of which a ne- 
gative resistance G, is connected. The long line is characterized by its wave resi- 
stance Z), propagation coefficient y = a+jf, and length /. It is assumed that the 
negative resistance has a nonlinear characteristic described by the equation 
(2.1) i= —S,u+S,u?, 


where 7 is the current and u — the voltage on the terminals of the negative resi- 
stance S, and Sj; are positive coefficients. Using the method of the third harmonic 


CG ee Fig. 1. The scheme of a _ generator 
t Dee with a long line excited by a negative 
Ss ( wil resistance G, 


[2], [3], it is assumed that the voltage on the terminals of the negative resistance 
has the form: 
(2.2) u = U,sinwt+ U,sin3wt+¢z). 

It is known that in this case the first harmonic of the current flowing across the 
negative resistance is expressed by the following approximate relation: 


Shick = = : Says 
(2.3) iy © G S,U3—S,U— = m;S,U? cos | sin wt— = m3S3U® sin y3 Cos wt, 


where 


(2.4) mM, = —— 


is the voltage third harmonic content. 
It can easily be seen that the amplitude of the current component in phase with 
the first harmonic of the voltage is . 


= é) Hs Ey | Ps 
(2.5) i a 4 Sg U*—S,U——- msS3U* Cos P2> 


and the amplitude of the orthogonal component is 


a 3 aoe Ne 
(2.6) vA — 4 m3S,U2 sin P3- 


In complex notation, we have: 
(2.7) I, = [e+ 1i/ tJel2) , 


In our considerations we shall use the notion of the admittance of the negative 
resistance, defined thus: 
I, 


(2.8) Ye = Gee 
Uy 
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Fig. 2 shows the phasor diagram for the first harmonic of the current and the 
voltage. The positive sign of the exponent j(z/2) is used for 3 <0. In this case, 
the negative (dynatron) resistance has a capacitive susceptance. For 3 > 0 the mi- 
nus sign is used and the susceptance is of the inductive type. 


Fig. 2. The phasor diagram for f.20) 
the first harmonics of voltage and s 
current 
Bet 


The components G, and B, of the admittance are expressed by the approximate 
equations: 


fa 

(2.9) iy U, a 753010 —m, cos P3)—S1; 
I’ 3 = ; 

(2.10) jaf es U; = — qimsS3Uj sin 3. 


3. The Admittance of the Long Line 


It is known that the input admittance of a long line short circuited at the end 
is determined by the equation 


__ ,, al cos B1+j sin Bl 
ep Bi Brg Bi-+jal sin Bl ’ 


where Z, — wave resistance, a — damping coefficient, 6 — propagation coefficient. 
If the damping is not too strong (al < 1), the following relations hold 


ale wi Re 
G2) a= *ruil-% Sale 


oa RG @ R? 
Go) aS oy Ee|1— Re + aC? * BotL? 
eR GZ, 
(3.4) See ehis 


where R, G, L, C are the constants of the line with values expressed in suitable 
units per unit length. 
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In order to avoid complicated expressions the following simplified equation 
will be used: 


by 
(3.5) Zo Vie 


(3.6) Bx oy LC. 

The components of the input impedance of the line 
OH) ZR 
are expressed by the equations: 

al 
e®) Bo eee pl+oPl? sin? Bl 
i sin 26/(1 —a?/?) 

(3.9) Ky =8Z, 


x ° cos? Bi-+ a2/? sin? Bl 
The components of the input admittance 


1 


(3.10) Y,=> 


= Gy tIBye 


are expressed by the equations 


Sp 1 al 
w Z, atl? cos? Bi-+ sin? Bl. ” 


(3.11) G 


1 sin 26/(1 —a?/?) 
(3.12) Bor St 7 a cost pT Raine OLE 


If a shunt fixed capacitance Cy is connected to the input terminals the resulting 
susceptance is 

1 sin 261(1 —a?/*) 
2Z, al® cos? Bl+sin? Bl. 


(3.13) B.S oc 


4. The Admittance Equation 
The steady state vibration parameters may be determined from the admittance 
equation: 
(4.1) Yip Ly = 0. 


This is a complex equation, that is it can be split up into two equations, the con- 
ductance equation 


(4.2) G,+G,, — 0, 
and the susceptance equation 
(4.3) B,+B, = 0. 
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5. The Conductance Equation. Determination of the Amplitude 


Substituting in (4.2) the expressions (2.9) and (3.11), we obtain the conductance 
equation in the form 
1 al 


3 oe 
Sl! Sa é = 
GD 4 sUi(1—m, cos ps)— S14 Zy al? cos? Bl-+sin? Bl us 


From this equation the vibration amplitude can be determined, expressed by 
the equation: 


Ss 1 al 
(5.2) C= *  Z__al* cos® Bl+sin® Bl” 


3 
a S3(1—mg Cos #3) 


Usually al< 1, bl» a/2 and m,cos gy; < 1. In this case, the vibration ampli- 
tude is expressed by the approximate formula: 


as Sial/Z, 
(5.3) U, 2 V sariget ¢ 


This is the same equation as in the case of a generator with a parallel resonance 
circuit having dynamic resistance R,, because the equation 

al 1 
(5.4) Sacey ae ae 


determines the dynamic resistance of a short circuited long line for the quarter 
wave resonance. 

Introducing the notion of the coefficient of the negative damping defined by the 
equation 


I 
(5.5) & = (s.-5,| Zy = S,Z)—al, 


5.6 ae Ee, 
(5.6) 0, a st 


6. The Susceptance Equation. Determination of the Frequency 


Substituting in the Eq. (4.3) the expressions (2.10) and (3.19), we obtain the 
susceptance equation in the form: 
1 sin 2 Bl(1—a?l?) 


Ki 
2Z, al cost pl+sintél + O° 


3 ane 
(6.1) <EpER m,5,U? sin @3 


312 S, Hahn 


Substituting suitable equations and rearranging as shown in the Appendix 1, 
this equation takes the normalized form: 


pein oy —ail?)—3yk (w" cos? 35 +sint3y>) 


4 : 
—— «* sin| arc tg 


al 
6.2 = — 
ee) y. a?l? a | ape 1 sin 3yz(1—a?/?) : 
: é Eee 1 
V (are cos? 3y5 + sin? 275} q’l? cos” 3y5 +sin? is 
sin 2y 5 (1—a2/2) 
2 + yk 5 ee) 
a*]? cos’ y +sin? Va 
with the following notations: 
623)5 v= —- . @y—quarter wave resonance frequency; 
0 
C, —the ratio of the fixed input capacitance to the total capacitance 
(6.4) k=— : 
Gl of the line; 
al —the coefficient of negative damping 


2 oo 
EM a ed a’l? cos? 6l+sin?Bl depending on the frequency. 


The Eq. (6.2) enables us to determine the vibration frequency y = w/q@,. In 
view of the complicated form, the graphical method is of practical use. It is con- 
venient to draw two curves on one diagram, one representing the dependance 
of the susceptance of the negative resistance on the frequency [the first, most 
complicated term of the Eq. (6.2)], (cf. Appendix 1, Eq. I. 17), the other — the 
dependance of the susceptance of the long line on the frequency, the fixed capa- 
citance C, being taken into consideration [the last two terms in the Eq. (6.2)]. 


7. Vibration Frequency for C, = 0 


Let us consider first the case in which the capacitance C, can be disregarded 
(Cy~0). In this case ¢ = é  [cf. Eq. (5.5)]. If now we reject a2/? as small in relation 
to 1, and if we consider the range where sin®y(z/2)>-a?/?cos?y(z/2), the Eq. (6.2) 
will take the following simplified form: 


beck e2 2/24 Ler eyes sin} arc (ee 
3 4 2) 2al Tt 
G21) aan =O OF 
cos?3y 5 +a?/? sin? 3y > 


Nonlinear theory of a generator 313 


The diagram of this relation for a certain value of the parameters ¢, and al is 
shown at Fig. 3. It is seen that the curve expressing the dependency of the sus- 
ceptance of the negative resistance on the frequency has the same form as for the 
frequency discriminator. The intersection with the curve ctg y(z/2) takes place 
at the point y = 1, that is the generator operates exactly with the quarter wave 
resonance frequency, independently of the value e. In the ideal case under con- 
sideration there are therefore no nonlinear influences acting on the vibration fre- 
quency. 


B B, 
-B 
: al=0.069 
06 42294 
a, 
Cl 
04 
02 
0 o<- 
-02 
-04 
-06 


Fig. 3. The susceptance diagram of the negative resistance B;, and the Susceptance 
of the line B, in function of the normalized frequency y = w/w). Note: By is taken 
with the opposite sign 


This property of frequency stabilization was foretold by J. GroszkowskI in an 
analytic way by means of the method of harmonics, [1]. 


8. Vibration Amplitude and Frequency, the Input Capacitance C, Being Taken into Account 


The Eq. (7.1) and the diagram of Fig. 3 concern the ideal case of Cy) = 0. In 
practice, however, the influence of the fixed capacitance Cy should be taken into 
consideration, therefore we are interested in the solution according to the Eq. (6.2). 
In addition to the quantity « and the parameter a/, there appears in this equation 
a parameter k = C,/C/ characteristic of the influence of the capacity C). 
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If nonlinear influences are rejected, the frequency is determined, according 
to the linear theory, by the relation: 


1 4 
(8.1) yk i was = 0. 


Fig. 4 shows a diagram of the dependency of the frequency on the parameter k. 
If, however, the nonlinear action of the third harmonic is taken into conside- 
ration, the vibration frequency is somewhat higher then that following from (8.1) 
and is determined by the Eq. (6.2). It should be added that with a frequency higher 
than the resonance frequency the admittance of the line is of the capacitive char- 


28 
26 
24 


2,2 


0 02 04 06 08 40 4, 


Gy 

Cl 

Fig. 4. The relation between the frequency and the ratio of the fixed 

input capacity Cy to the total capacity of the line Cl; y; = w/a, is 

the normalized frequency of the fundamental mode, y, = w,3/a@ —the 
normalized frequency of the second order mode 


2 
k= 


acter and the admittance of the dynatron is of the inductive character. We can 
speak of the «inductance» of a dynatron in place of the notion of the «<capaci- 
tance» of a dynatron introduced by GROSZKOWSKI in the case where a circuit 
with fixed constants, [2], is excited by a dynatron. The physical cause of the dif- 
ference lies in the fact that a long line has an inductive admittance for the third 
harmonic, while the circuit of parallel resonance has a capacitive one. 

Fig. 5 shows graphical solutions of the Eq. (6.2) for various values of the pa- 
rameter k. It is shown that with increasing k, the susceptance curves of the negative 
resistance are shifted to the left in relation to the curve for k = 0. In addition 
they are asymmetric and show a decreasing slope in the middle portion. It is 
essential, that the intersection of the curves for the two susceptances takes place 
below the axis of abscisae, which means that the susceptance of the negative resis- 
tance is of the inductive type and that of the line of the capacitive type as explained 
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above. Since the scale of the curve for the susceptance of the line varies with varying 
9, the frequency varies somewhat with varying é). This variation depends on the 
slope of the discriminator curve in its middle portion. It is seen from the Fig. 5 
that this slope decreases for large k. However, in view of the great slope of the 
curve, the frequency variation due to the variation of €) is much smaller than 
in a dynatron generator with a circuit with fixed constants. 


B, i~Byh 


x eg Oss i ol =0.069 
moNe ae 
ESR Os \ 3 
fhe ee ss | Sere 
\ \ 1\ E 

oe 1 \ ——-— k= 
oat- \ \ : \ bbs tretses k=03 
aX ty eee 08 
ce : Siena hed D 


~. 


yf SSS 


a N poe Ss 14 12 x re a 
-02 x . l : 
s oe 
04 s : 
- 06 
—08 


Fig. 5. The diagram of the susceptance of the negative conductivity B, and the sus- 
ceptance of the line B, in function of the normalized frequency y = w/w, for various 
k = C,/C;. Note: By is taken with the opposite sing. 


It is characteristic that above a certain limit value &k the curves of Fig. 4 do 
not intersect, which means that the phase condition cannot be satisfied although 
the negative damping coefficient « is greater than zero. For, according to the Eqs. 
(5.2) and (6.5), this takes place when 


al 
(8.2) PSL) , = \e50: 


a?l? cos” IS +sin? y> 


If sin?y(2/2)s>a7l?cos*y(a/2) we can write in an approximate manner 


(8.3) e= |S 7. ae 0. 


6 Problemy drgan 
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This equation can be rewritten thus 


(8.4) i € = &+al po Ae >0. 


sah, 
sin’ y5 


©The relation ¢ = f(y) is drawn as an example in Fig. 6 for al = 0.069. The bounda- 
ry condition «=0 corresponds to the value of frequency determined by the 
equation: Sie 


(8.5) sin py Se 


Fig. 6. The dependency of the negative damping 
coefficient « on the relative frequency y = w/w. 
The damping al = 0.069, s,z, = 0.343, ©,=0.075. 


It is seen that in the concrete case represented by Fig. 6 the amplitude con- 
dition is satisfied for y>0.3, while the phase condition is satisfied if y>0.75 (ap- 
proximately). 


9. Generation of Vibration of Higher Orders 


A generator with a long line has theoretically an infinite number of degrees. 
of freedom. In the ideal case of a line of zero input capacitance and with the dam- 
ping independent of the frequency, the generator can, according to the linear theory, 
excite itself on every frequency, for which 
(9.1) ctg yo = 0; 

Z 
that is, for y= 1, 3, 5,.... The steady state frequency of the generator would be 
accidental. If the influence of the input capacitance is taken into consideration, 
the frequencies of vibrations of higher order are not in a harmonic relation and 
are determined by solving the Eq. (8.1). For the nonlinear effects the Eq. (6.2) 
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a) BBW b 


b) 8 -By 


B 


By 


i 2 \| 3 _w 
Wo 


<4 


Fig. 7. Diagram of the susceptance of the negative resistance By; and 
the input susceptance of the line B, in function of the normalized frequency 
y =@/wo in the range of the fundamental and the second order mode: 
naa — 5); — 01250 (D) PA 10:3, 
6* 
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can be used by enlarging the frequency range in an appropriate manner, for which 
the susceptance curves are determined. Fig. 7 shows the curves for a second order 
vibration for the parameters k = 0.2 and k = 0.3. Itis seen that in the given 
case the intersection of the curves takes place at three points. The problem as 
to whether any of these points determines a stable of vibration remains open, 
however. This problem will not be considered here. We shall only point out the 
fact that the stability condition has the form, [3], 
0G OB OG OB 


OU Gp aay OU 


where 
'G=G,+G, B=B,+B,: 
[see Eqs. (4.1), (4.2), (4.3)]. 
10. Conclusions 

On the basis of considerations taking nonlinear effects into consideration, the 
mechanism of steady-state vibration in a generator with a long line short-circuited 
at the end is described. In particular, it is shown that in the case under conside- 
ration the dynatron has an inductive admittance. In addition, it follows from the 
considerations that in the case where the input capacitance of the line is too high, 
the phase condition is not satisfied in spite of a negative damping of the system. 
Unfortunately, this fact has not yet been confirmed experimentally. The chief dif- 
ficulty consists in that the ideal case of negative resistance of cubic characteristic 
was considered. The characteristic of the nonlinear resistance used in the experi- 
ments was asymmetric, that is it can be described by an equation containing a 
quadratic term. In consequence, there was an influence of the second harmonic. 
The investigations of the experimental layout being of a qualitative character, 
neither a description of the layout nor the results will be given. However, the 
qualitative investigations confirmed the conclusion that due to the influence of 
the third harmonic, the generator produces vibrations of frequency higher than 
the resonance frequency of the line. The influence of the second harmonic causes. 
a decrease in the frequency of the generator, the line admittance of the second 
harmonic being of the capacitance type. 

Generators with long lines are most commonly used in the uhf aa the micro- 
wave range. Their electrical layouts are more complicated than that in Fig. 1. In 
addition the negative resistance characteristic may have, as was mentioned, a shape: 
different from that used in the considerations. For these causes, the theory given 
here concerns a certain ideal case and its object is a qualitative explanation of 
nonlinear influences in generators with long lines. Let us observe, however, that 
the equivalent scheme layout of a generator with a long line excited by a tunnel 
diode approaches the scheme in Fig. 1. The characteristic of the negative resistance 
of a tunnel diode is of the dynatron type and may be approximated by the charac- 
teristic according to the Eq. (2.1). The extension of the present theory to the case 
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of a characteristic containing the square term is possible but the equations become 
very complicated. 

A generator with a long line has many degrees of freedom and investigation 
of the possibility of self excitation of higher order vibrations is also essential. In 
principle, the Eq. (6.2) enables the investigation of the characteristic of ‘the steady 
State also for higher order vibrations. In practice, however, the fact that the dam- 
ping of the line is a function of frequency should be taken into consideration- 
In the experimental generator, no higher order vibration was observed either 
in the steady state or during build up of vibrations. 

The theory presented was the subject of a lecture held by the author at the 
conference on non-linear vibrations organised by Czechoslovak and Polish Academy 
of Sciences at Liblice near Prague on September, 21-24, 1960. 


Appendix 


In agreement with (2.10) the susceptance of the negative resistance is expressed 
the equation: 


3 eer 
(A. 1) B, => aMsS3U" sin P3- 


Since the dependency of B, on frequency is desired, therefore we must know 
in which way m, and @, vary with the frequency. 
The current of the third harmonic is expressed by the approximate formula 


(A; 2) i 5,0" sin wt. 


Its amplitude is therefore equal to: 


= (Bpeee 
(A. 3) : Ts = 7 S3U*. 


The content of the third harmonic is determined by the equation 
U. 3 J 3 Zsy | 


(A. 4) Ms; = 0, = Ui, 


where |Z;,,| is the modulus of the input resistance of the line for the frequency 
of the third harmonic, the capacitance C, connected in parallel, being taken into 
account. 
The input conductance of the line is expressed by the Eq. (3.11): 
is al 
eo Cow = Z, al? cos? 381+ sin? 361° 


The resultant input conductance of the line, taking the capacitance C, into 


consideration, is 
1 sin 6pl(1 —a?]?) 
2Z, al? cos? 361+ sin? 361° 


{A. 6) e306, 
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Hence the modulus of the input impedance of the line is 
1 
A, 7 Zé, Ww — SSS SS SS 
1 


- 1 oP? ei er es 
Z? (a? cos? 3f/-+ sin? 361)? +| ONO 9Z,, al? cos? 3B] + sin? 361 
Substituting (A. 7) and (A. 3) in (A. 4), we obtain: 
| 5 Si02Z 
(Ais 8) tgp ee ee 
al? + [30c,z 1 sin6fli—ol) - 
(el? cos? 3B/-+ sin?361)? +| OX 00 Fel cos? 3B1-+ sin? 381 


The angle v3 will be determined by taking the relation 


(A. 9) Us = —J3Z ay; 
where 
(A. 10) LZ sy a R3y ads au 


is the impedance for the third harmonic. Let us construct the phasor diagram 
for voltages and currents at the moment ¢ = 0 (Fig. 8). 


Uy 
inductive(x,>0) capacite (x;<0} 


Fig. 8. The phasor diagram of voltages and currents of the 
fundamental and the third harmonic 


It is seen from the diagram that if X,;>0 (the line shows the inductive pro- 
perties for the frequency 3w), than the angle @ 3 is positive and contained in the 
range 
(A. 11) 1 /2Z=P320. 


If X;<0 (the line shows the capacitive properties), then the angle gv is negative 
and contained in the range 


(A, 12) —1/2<y3<0. 
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The angle v3 will be found from the equation 
38w . 


A. 13 t eee ee OW. 
( ) £ P3 Gs ) 


or, substituting (1.6) and (1.5), we obtain: 
1 sin 66/(1 — al?) 
IZ, aPcos!3 I+ sint3 1 °C 
is al 
Z, Il cos? 61+ sin? 61 


After some transformations, we obtain: 
ie 
— sin 661(1 —a?/*) — 3wC,Z,(a*l? cos? 361+ sin? 361) 


(A. 15) Y3 = arc tg 
al 


Substituting (A. 8) and (A. 15) in (A. 1), we obtain: 
1 
— sin 6f1(1 —a?l?) — 3wC,Z,(a2l? cos? 361+ sin? a 


(A. 16) 
_3 s2sz, sin [are tg 
pee 4 al 
= 
ol? 1 sinopld—a@r) 7 
V om cost 3A1-+ sin? 381% * [2c.20— 2 “oPcos*3Bl-+sin236l | 


Substituting in this equation the relation (5.6), we obtain: 
mt sin 661(1 —a®l?) — 3mC,Z,(a2/? cos? 36] + sin?3Bl) 


(A. 17) 
1 
Steg? soi 
mee 3 © Zo Sin arc tg A 
a ere + 30c,z 1 sin 661(1 — al?) a 
(a?/? cos? 361 + sin? 361)? On 00D GI cos? 361-+ sin? 361 


The susceptance equation (4.3) and (6.1) takes, on substituting (A. 17), the 


form: 


(A. 18) _ sin 61(1—a2I2) — 30C,Z, (a2? cos? 361+ sin? 361 
ae PIT 
3 é* sin | arc tg ay “j 
oer | sn6pld—e@y) 2 
(a2 cos® 380-4 sin® 3p) | [2CoZo— 2 oP cos? 3p1-+ sin® 3f1 | 


1 sin 261 (1 —a7/?) oe 
2 al? cos? BI -+ sin? fl ee aaa 


This equation may be reduced to the normalized form by denoting: 


> 


ol oe ay where y= es ee 
wager 5? Wo ma eV 
Co 


@ M4 
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We have assumed Z,=J/L/C. With these notations, we have: 
(A. 19) 


uo! 


a3 
é“ sin j arc t 
8 al 


a*l? 1 sin 3yz (1—a?l*) 
a ae ee 


Ol Cas ay + sin? 3y 


5sin ahah! — 3yk (wr cos? ay + sin? 3y 5 al 
7 
> 


V (or cos” sy + sin? 375] 


1 sin yx (1—a?/?) 


GC isha a 
a*/? cos" y-~ -— sin’ y-> 


For the derivation of the Eqs. (A.17) to (A.19) it was assumed that the amplitude 
of vibration is expressed by the Eg. (5.6) and does not depend on fi. If (5.6) is 
replaced by (5.2) then in the Eq. (A. 19) the following expression should be substi- 
tuted in place of the constant coefficient of negative damping 


2=(S,Z,— = : 
~ \ Fa]? cos? BI+ sin? Bl] ° 
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Strieszezenie 


TEORIA NIELINIOWA GENERATORA O LINII DLUGIEJ 
POBUDZANEJ OPOREM UJEMNYM 


Rozpatrzono generator zlozony z lini dtugiej Cwier¢falowej na koricu zwartej, 
pobudzanej do drgan oporem ujemnym dynatronowym o charakterystyce szeé- 
ciennej. 

Stosujac metode rownania admitancji utozono réwnania, z ktérych mozna wy- 
znaczy¢ amplitude i czestotliwosé drgan. Rozpatrzono przypadek idealny z pominie- 
ciem wplywu pojemnosci skupionej C,, podtaczonej rdéwnolegle do zaciskéw 
wyjsciowych linii, jak réwniez przypadek z uwzglednieniem wpltywu tej pojem- 
nosci. 

Wykazano, ze jesli pojemnosé jest zbyt duza, to warunek fazy nie jest spetniony, 
mimo ze odttumienie generatora « jest dodatnie. 
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Pesrwme 


HEJIMHEMHAA TEOPHA TEHEPATOPA C JJIMHHOM JIMHMEN 
BOSBYXKTAEMOM OTPHIATEJIBHbIM COMPOTMBJIEHHEM 


PaccmaTpuBaetca reHepaTop, COCTOAMIMM U3 KOPOTKO-38MKHYTOH UeTBEpTb-BON- 
HOBOH JIMHMM, BOSOyKTAaeMOM OTPHWaTeUbHbIM COMPOTHBIICHHeM WMHATPOHHOTO 
Tuna, C KyOwuecKOH xXapakTepuicTuKoN. 


IIpumeuaa MeTO ypaBHeHHA aMUTaHCa COCTaBJIAIOTCA ypaBHeHHA, 10 KOTOPBIM 
MO?KHO Ope CIMT aMMIMTyAy MU uacToTy KoeOaHui. 

PaccmaTpuBacTcA WpeanbHbI ciyualt upuH mpeHeOpexeHuu BIMAHHeM cocpe- 
WOTOUCHHOM CMKOCTH Cy, NOAKIOUCHHOM NapasIeNbHO K BXOJHbIM 3a>9KUMaM JIMHHN, 
a TakoKe cylyuelH c yudeTOM BIIMAHMA 9TOH eMKOCTH. 

Jloka3biBpacTCA, UTO CCIM €MKOCTb CMUIKOM OobUIad, Tora ycroBue (a3sbI 
He YOBJICTBOpeHO HeCMOTpA Ha TO, UYTO OTPHMaTeNbHOe 3aTyXaHHe KOeOaHuit 
reHepaTopa € ABIIACTCA MOJIOYKUTCIIBHBIM. 
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THE UNLOADING WAVE IN A LAYERED BODY 
WITH RIGID UNLOADING CHARACTERISTIC 


SYLWESTER KALISKI (WARSAW) 
1. Introduction 


The problem of propagation and reflection of a loading and unloading wave 
from a rigid and elastic wall was solved, for a body assumed to have a rigid unlo- 
ading characteristic, in Refs. [1] and [2]. In spite of the qualitative difference 
between this solution and the reality (in the unloading zone, the body behaves 
as a rigid body subjected to an initial deformation) the quantitative results very 
closely approach the reality. A similar problem for a free surface was treated 
in the paper [3]. Analogous solutions for non-homogeneous bodies with non-rigid 
unloading characteristic were investigated by using the technique of numerical 
solutions in J. Osteckt’s doctorate thesis, [4], where the limit cases of solution 
for rigid unloading characteristics were discussed and where the above statements 
found confirmation. 

As regards practical problems, we are essentially interested, however, in solutions 
convenient for a designer — in a closed form, if possible. In favour of this interest 
there is also the fact that approximate mathematical operations result also in an 
error, which generally increases with time if the method of nets is used, (cf. [4]), while 
for a physical approximation, for which an accurate mathematical solution can 
be obtained, the causes and the extent of error are fully known. For these reasons 
it appears in the author’s opinion, to be worth while, to consider some other so- 
lutions, based on the assumption that the body has a rigid unloading characteristic. 
Since the problem of plastic waves in layered bodies is often met with in practice, 
a solution for such a body will be presented. 

The problem of the unloading wave will be treated, that of the loading, wave 
being well known and considered to be classical. 

In subsequent papers, the problem of the unloading wave will also be solved 
for a cylindrical wave and a spherical wave — which requires another approach, 
in.view of the composite state of stress. 
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The solution given below can easily be generalized to the case of bodies non- 
homogeneous in each particular layer and to cases of curvilinear or angular loading 
characteristic (Fig. 1). 


ron os 
0 
of) 
Os 
E 


Fig. 1 Fig. 2 


However, for the sake of clarity, and in order to obtain relatively simple equations, 
we shall consider the cases of the (o—«)-characteristic represented in Figs. 2 and 
3. Both cases find practical application to various bodies, in soil mechanics in 
particular. For moderate soil pressures, diagram 3 may sometimes be assumed 
as an approximation of the real diagram given in Fig. 4. 


(o} 


Fig. 3 Fig. 4 


The viscosity is not taken into consideration because we shall confine ourselves 
to pressures at the soil surface varying relatively slowly in time in the unloading 
range (certain explosions, for instance, for which the unloading time is of the order 
(0.1—1.0)/sec). It is of course an approximate assumption. 


In addition, in order not to complicate the solution for reflected waves, we 
shall consider, similarly to [2], a strong discontinuity wave. Weak discontinuity 
waves involve, for reflection, certain additional complications which do not, howe- 
ver, change the essential features of the problem, the solution thus being con- 
structed following the rules given below. Therefore this problem will not be treated 
in detail. 

Let us consider first the general solution in a qualitative manner, by means 
of the phase plane, for the (o—e)-characteristics as shown in Figs. 2 and 3. 


Then, some simple examples will be given. As already indicated, we can con- 
struct, in principle, the solution for a (o—«)-characteristic as represented in Fig. 1. 


———— ELL 
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The solution becomes considerably complicated, however, and loses clarity be- 
cause of, among others things, the fact that the unloading wave becomes a curve, 
[1], even for waves of strong discontinuity. In view of the requirement of simpli- 
city this case will be omitted. 


2. The Unloading Wave in the Case 0.0, > 0,0, 


Let us consider the one-dimensional case of the semi-space to which pressure 
is applied suddenly on the surface, and then decreases monotonically. A two- 
layer medium is assumed to occupy the semi-space. In the case of a greater number 
of layers, the solutions will be repeated. (Then, there exists also a possibility of 
the appearance of weak discontinuity waves). 

Let us consider three cases: (1) 6,5 =0,.=0, (2) |o,|<|o49|, (3) |oa| >|o,9|. Although 
the first case is a particular case of the other two, the relations in the phase plane 
change in this case in an essential manner, therefore it will be treated separately. 
Let us start from this case. The (o—¢)-diagram is that of Fig. 2. 

Case 1. 6, = 0,2. = 0. Since various sub-cases may be considered for various 
durations of the pressure on the surface, various distances of the second layer etc., 
therefore only one of them will be considered, other variants which can be solved 


P(t) F, 


Fig. 5 


in the same manner being disregarded. The parameters will be assumed as shown in 
Fig. 5. Each particular solution zone is denoted by a different number in Fig. 5. 


The relevant solutions will be denoted by the same labels. 
In order to fix the attention, a, >a, has been assumed in Fig. 5, which corres- 


ponds to the condition Qa, > 04, for 0, = Q2. 


7 Problemy drgan 
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Let us establish the equations and obtain the solution for each particular zone. 
We assume that P(t) is positive for compression. 
Zone I. In agreement with the results of Ref. [1], we have: 


1 
01Qt 


a x 
(2.2) =— Ex H | Fin 


(2.3) a= : [ea ea a Jen 
xy ay xy ay 


(2.1) i 


H(t) Ear Vip» 


3 a 
oe ay 
(2.4) O1p a 2 ) 
where 
(2.5) HQ) = { P(t)dty, 
0 


and the subscript p denotes the values of the functions at the front of the wave. 
Zone II, The equation of motion in the Zone II is 


dy 
(2.6) 0,(2/—a,t oy — P(t)+04), 
with the initial condition 
1 
P| Vv. = ——-A(t,) = v1 (t for f=; 
(2.7) 1 Sait) = ne) ; 
In addition, 
(2.8) Eg = E1p. 
The stress o,(x, ¢) can be expressed thus 
(2.9) d(x, t) = 0,,(x)— Ao,(x, 0), 


where the symbol Ao is explained by Fig. 3. 
For an element dx of a bar AB, in the Zone II of the phase plane (Fig. 5), we 
can write the equilibrium equation 


Ove O0_ 
(2.10) Our ae 0, 
or 
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with the boundary condition 
(2.12) Ao, = y(t) 100 Xi == 05 


where y(t) = P(t)—P,. 
Hence 


x 


¥ 0064p OVe 
(2.13) Ao» ji wane st) t p(t). 
0 


Thus all the quantities of interest have been determined for Zone II. 
Zone V and VIII. \n the zones V and VIII the velocities of bars of unit cross- 
sectional area separated mentally from the semi-space should be equal; therefore 


(2.14) eae. 

The stress excess o, at the front of the reflected wave is 
CAE Op = —014,(V2—Vs) 

and 

(2.16) Oxy = Og = 01,0. 


The equations of motion determining vy; and y, [for v, determined from the Eq. 
(2.6)] have the form: 


dy, 
01(a,t—l) a = —(09+01,) +6, 
(25137) As 
Q2[a2(¢—t,)] = —O,+Ogr,= —O;— 0242s, 


or 


dv 
(2.18) [(0141 + G22) t — (01! + 02427) ] ae + (0141 + 0242)V5 = — yp + O12, 


with the initial condition 


1% 

: 7 ——_—. 
(2.19) V5 (73) Vy (Ty) 014, + 004, 

For 0, = @2 =e the Eqs. (2.18) and (2.19) take the form: 

d 
(2.20) ol (ay+as)t—(I-+-asr)) 52 +-0(ar + 42)¥5 = Fp + OarYe > 
a 

(2.21) v,(T3) —e, Ze) rae . 


Next, we shall have of course: 


(2:22) a are 
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o; and o, are computed in the same manner as in the zone II, except that for 
the determination of the function y(t) we have now the conditions: 


(2 23) 05(x, t) = O1p(X; t)+0o(x)— Ao, (x, t), 
Kop = 0 for x= 2/l—a,t; 

(2 24) 04(x, t) = Ogy(x)— Ao, (x, t), 
One 0 for x =I1+a,(t—t). 


Zone III, VI, and IX. We have 
(2:25) Veer 


The equation of motion is 


dy. 
(2.26) [erail+Q2a2(t—11) > =P(t)+05,, 
with the initial condition 
(2.27) ES Tp for hee 
(2.28) Es = Ep, Ee = Esp Eg = Egy: 


For the remaining zone, o,; and o will be computed according to formulae ana- 


logous to those for Zone II, except that for the determination of y(t) we use the 
condition 


(2.29) Ao; = Ao, = 0 for x= 2/—a,t, 
Noa for x =//+a,(t—t). 
Zone IV, VII and X. We have, as before, 
(2.30) Va Va — ies 
dy, 
(2.31) [o14/+ 0242(t—T) oe O10: 
with the initial condition 
232) TAC) Sia 
Next, 
(2233) &4 = €1)» &7 = Esp, E10 = Egps 


o will be computed as before, with analogous conditions for Ac, from which 
y(t) is calculated. If we had t, >t or t* > 1, the configuration of the zones would 
be different; the principle of constructing these solutions would, however, remain 
unchanged. 


Let us consider now the second case, that is 

Case 2. logl < lo, 1. Similarly to the first case various solution variants and 
various zone configurations will be obtained depending on the assumed parame- 
ters T, 7, 1, 043, O,. and P,. We shall not, in view of limitations of space, discuss 
all the possible variants. We shall confine ourselves to a single case of essential 
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importance. The principle of constructing the solution will not be changed for 
another set of parameters, only different qualitative and quantitative conclusions 
will be obtained. To fix the attention, the set of parameters will be assumed as 
shown in Fig. 6. 


P(t) P, 9s, Xo l 


Fig. 6 


The following solutions are obtained for each particular zone. 
Zone I, The Eqs. (2.1)-(2.5) preserve their form except for H(t) now replaced by 


(2.34) H() = (P(t) +o,,]dh, 
0 


and for the computation of o the value o,, should be added to the right-hand mem- 
bers of the Eqs. (2.3) and (2.4). 
Zone II. The equation of motion in Zone II is 


dy 
(2735) 0,(2/ ayt) or = P(t) | O1p> 


with the initial condition: 


(2.36) Vo(T) = Oat, H(t) = v4(%). 
Next, 
(2.37) &5 = Erp» 


where, in view of a different definition of H(t), €, and o,, have sense other than 


in the case o,, = 0,,=—0. 
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To calculate o,(x, f), equations are available analogous to those for the previous 
case 


(2.38) 02(x, t) = 0;,(x)—Ao,(x, 2), 
* 1 0 oy. 
(2.39) Nop = i | ee 01 ai} | P(t)—Py. 
0 
Zone V and VIII. We have 
(2.40) V5 = Vg, 
(2.41) Op = —014,(V2—Vz). 


The equations of motion are: 


d 
ox(at—) 8 = —(09+01,) +0, = 0, +010 (V2— V5) —A 1p, 
(2.42) 


dy 
02[a2(t—T,)] 7 = —O, +O gp = 01 —0202Vg + O52, O52 <0, 
or, eliminating o, and making use of (2.40) 


(2.43) [outs + 0003)t— (oul + over ) 7+ (01d1+0202)¥5 = = —(01,—O 52) +0141V2, 
with the initial condition 


(2.44) ¥;(t3) = 2v,(7) —2 Ose Ost 


0144+ 0242 0141 +0242 ; 


In addition 


(2.45) & = ea ‘ —— 

and 

(2.46) Bee t)= 01)(x)+09(x)— Ao, (x, t), 
d3(x, t) = o,(x)— Ao, (x, 2), 


with the conditions for the determination of y(t), analogous to (2.23) and (2.24). 
Zone III, VI, IX. In the case of t > ¢°, where t° — the time for which o,, = o,9, 
we have: 


(2.47) Ve = vp = O, 

(2.48) &g= fp, && = 5p, € =f, (& =0 for. "7s 1,)- 

v3 will be calculated from an equation identical with (2.35), and the stresses accord- 
ing to equations analogous to those for Zone II. 


Zones IV, VII and X. For t >t*, where ¢* is the time for which o, = 0, we 
have 


(2.49) V4 = V7 = Vio = O, 
(2.50) y= 1p, & = Exp, Erg = Esp, [ey = O {Of Wes, ], 
(291) 0, = O7 = Oy = —P(t). 
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For t >t, we have o = vy = 0, and permanent deformations take place as in Zones 
IV, VII and X. 


Case 3. |o,|> |o,.|. The parameters are assumed as shown in Fig. 7. 


———— 
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P(t) P, Os, i l ~ 


Fig. 7 


In Fig. 7, it has been assumed that the pressure increases in the interval /, and de- 
creases starting from t,. The case where for t) < t, the pressure reaches its maxi- 
mum at the point / can occur only if P(t) decreases in such a way that the plastic wave 
in Zone II ends before x = / is reached. Then, in the medium, a wave of weak dis- 
continuity appears (Fig. 8) [1]. 


i 


Fig. 8 


This case is somewhat trivial and will not be discussed here. The possibility of ap- 
pearance of such a maximum for t% > 7, will not be discussed either, in view of 
the confined scope of the present paper. 
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Let us describe the case of Fig. 7, where 6,,,, takes place for t,. In the discus- 
sion of the equations for each particular zone, the equations will be written for 
v only. For, having v, we are able to determine o and « according to the above 
prescriptions. 

Zone I. The equation of motion has the form: 


dy 
(2.52) e(l—ayt) = 0 +0; = —O4 FO s2— 02021, 
with the initial condition 
(253) y,(0) = 0. 
Zone IX. 
(2.54) DES, ( a =| ' 
a 
Zone II. The equation of motion has the form; 
dy 
(2.55) Quaat = P(t)+ on, = Pt)+04—01%(V2— V7), 
with the initial condition 
(2.56) (0) aa 
Q1% 


Zone III. The equation of motion has the form: 


dy 
(2.57) 0,(21—a,?) ar = P(t)+o,(2/—a,t) . 
with the initial condition 
(2.58) ¥3(T1) = V2(T}). 
Zone VI and X. In these zones, we have: 
(2259) Ve = Vio» 
(2.60) Oy = —014;(V3—Ve), 


dy 
ex(at—l)— = —(6)+0,,)+0, = 01+ 0141 (V3— Vg) — O25, 
(2.61) 
Ao 
0242(t—T}) Wik = —O, +19) = —9,+-G,2—0242(Vip— Vo), 


or, eliminating 0), 


d 
(2.62) [(0141 + @242)t— (0,!-+ 09457 )] E+ (ers +0002)¥ = 


oe — (G2p—O 59) +0141 V3+0242V9, 
with the initial condition 
(2.63) Vg(T ) —— V2(T1) Q1% = O52 —O4 + 0191V1 + 02A2Vq 


019, + 0242 0141+ 0242 


> 
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= 


where 
Vo =, (7). 


In the remaining zones we obtain, bearing in mind the changes in numeration 
of the zones, the same relations as in case 2. 

In view of the limited scope of this paper, the discussion of other variants for 
various values of the parameters Po, 6,4, 0,2, / etc. assumed will be omitted. 

In spite of different qualitative and quantitative results for different values of 
these parameters, the general approach will remain unchanged. 

The choice of the parameters for the above cases has been so made that the most 
essential cases are discussed, the ones which are trivial or similar to those discussed 
above being omitted. 


3. The Unloading Wave in the Case of 0,0, > 0.4, 


Three cases should be considered, as in Sec. 2. 


(1) On = O52 nae 0, (2) O51 <= O52], (3) lo a |o,2]- 


Here also a considerable number of variants is obtained, which cannot be discussed 
within the framework of a brief paper. However, since the essential features of 


th 


ae = 
x 


Pit) P, 0 l 
Fig. 9 


the method remain unchanged, as compared with Sec. 2, the discussion of the cases 
(2) and (3) will be omitted, case (1) being the only case discussed as an example 


@ase lid. — 0.5 — 0. 


Let us assume that the parameters in the phase plane are as shown in Fig. 9 
where a, < a, has been assumed, which corresponds to the condition 


014, > 0242 for 01 = 02 = 2. 
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This case is much simpler than the analogous case of Sec. 2. Let us consider now 
the equations and the solutions for each particular zone. 

Zone I. The equations of Sec. 2, for the case of o,, = 0,. = 0, preserve their 
validity. 

Zone II. The equation of motion for Zone II is the equation for the unloading 
zone, and the stress at the point O’ will undergo a jump-like decrease. This equation 
has the form: 


d 
G.1) [ol-+ova(t—ty)] 2 = P()—0uaaves 


with the initial condition: 
(3.2) v2(T,) = vz (7). 
The distributions o, and ¢, will be calculated similarly to Sec. 2. 
Zone III. The equation of motion has the form 
dy. 
(3.3) [01/-+0242(t—T})] ar + 02G2V2 = 0, 


with the initial condition 
(3.4) v3(T) = V2(T). 


o3 and ¢; will be computed in the same manner as in the foregoing section. 


4. Example 


Two examples will be considered. In view of the great number and variety of 
cases, and the limitet extent of this paper, let us consider only some examples for 
the case 1, that is for o,, = o,. of the Sec. 2 and 3, in view of their simplicity and 
practical importance in soil mechanics where the simplified diagrams can be 
applied to approximate the real (o—e)-diagram (Fig. 4). 


Example 1. 
Let us make the following assumptions: 
01 = 02 = @, ay < a, 
t 
(4.1) P(t) = Ps(1 al for BP ea 
Plt 0 for. #7, 


Then, in agreement with the considerations of Sec. 2, the solutions for each parti- 
cular zone are as follows. 
Zone I. 


(4.2) ie A= (! seat 
1 


(4.3) eee ee = op 
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4.4 i Pigalle 

P, ay P, a 
4.5 i o>] e = 
>) eas EF; | Te) oa? (! wat | 


Zone II. The equation of motion (2.6) becomes 
dy t 21—ay,t 3) as 
4.6) ola)? = (1 —7}-Al! ae 7 —P(5 = 


with the initial condition, (2.7): 


IPs f 
(4.7) 2(t) = 4(%) = mAb = | ) 
ea, De 
On integrating, we obtain: 
IS ae) eg t Set 1 T 
4.8 = 2} ——}= iP 
8) ih a fe (2 Sle SANG | 
Next, 
iP, x 
4.9 = ——2(1-— 
Ce) 2 oar (! fe, 
and, on the basis of the Eqs. (2.9) — (2.13): 
x t 2x t—T 
=—P zl Na 
(4.10) HO ar cre oreary aeaed 


Zones V and VIII. We have v; = vg and the equation of motion for both zones 
(2.20) is 


d 
411) ol(ay+a.)t— C+ ayt,) |? +0(@,+.44)5 = — oy +a," = 


i 2/—a,t T eat Seles ott |e 


2a,t 1 


= 2P,| 71n(2 ange he : a 
i} 


\ 


or, after transformation, 


dv; seers T ase t Eviaan WB 
(4.12)  (t—t) peaet: ene one ean 5 2 |=0. 


The solution of the Eq. (4.12) may be represented in the form: 


(4.13) iO a J Maar, +C| 


From the initial condition (2.21), it follows that C=0 because for t— 1}, 
in the expression (4.13), we obtain, for C = 0, 


(4.14) y,(ty) = eho es: 


o(a, +42) 
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which coincides with the initial condition (2.21). Writing out (4.13) we obtain 
after some manipulation: 


ZF, ed 
0 (a, +42)(t—T}) | 


Sats Ep Tet 
+(I 5 aK T1) maces t)In 2 : 


(4.15) va(t) = v5(t) = 


Next, 
1 21—a,t ty CC Se 
419 da ln 2 =O anne ye et 
Nes {2 5 Tt T% a 
t —— (27,—A)l : 
0(a,+a,)(t—T) pads aK T) t ( Gq In % at 
—] 
(4.17) = 


a, 
Similarly, by virtue of the equations of the type (2.9)-(2.13) and making use 
of (2.23) and (2.24), we obtain: 


[j= Dad ot i mal 
(4.18) v0) =—Poll z p+a(2i sae ri eae 


2ayt i aT, T 

2a, e—t? ow oan iA 21—a,t cs 

Centerline i: ay ee 
tse (ato 
1 2 
oO 2P, ieee > Shula 
(4.19) a,(x, t) (CE ee i| aE ie Ds ee T) 
t 


2 Qr, in | at | ire es aes i-20G= “ot. 


For a, = oo, that is for the case of incidence and reflection of the unloading 
wave in the first layer from a perfectly rigid wall, we obtain 


Oe: t 6 ale Py 9G Ah t 
(4.20) = 274] 1+ t a4 in(2—"). 


2 1 


This formula coincides with the formula obtained in Ref. [2] for the discussion 
of the problem of reflection of the unloading wave from a rigid wall. 

Zones III, VI and IX. In what follows, we shall give only the computation of 
y. Having v, the quantities « and o can be computed as before. Since v3 = Vg = Vo, 
in agreement with (2.25), therefore the equation of motion (2.26) has the form 


: 
(4.21) olla Ge = Pll -) Heat 
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with the initial condition 
(4.22) v3(t*) = v2(t*). 


Integrating this, and making use of the initial condition (4.22), we obtain, after 
rearrangement: 


P, t Dele it 
4.23 tt : 
42) WO agony z)ta[ eB -el+ 
CE ais r ehaets 
pty ec et (1-5 
The value of ¢* will be computed from: 
(4.24) Vo —— V5. 
It can be computed by assuming a numerical value of 7,/z (in view of the character 


of the Eq. (4.24), usually transcendental). It has been assumed that 1, < t* <t 
and t, >1/2, in agreement with Fig. 5. 


Zones IV, VII and X. In agreement with (2.30) and (2.31), we have 


Va a7 
and 
dy, 
(4.25) ell-+a,(t—t)) +0a,v, = 0, 
with the initial condition 
(4.26) V(t) = ¥5(7)- 


The Eq. (4.25) is identical with (4.21) if it is assumed that P(t) = 0. Integrating, 
this, we obtain: 


(4.27) v4 (t) = P a shes E n2— T+ 


U+a(t—ty)lo (2 a 


= 


Ye . 
eee at iiner| 


The above remark on ¢* remains valid. 
Thus the full solution is obtained in a closed form for the case considered. 


Example 2. The data are the same as in the Example 1, except that a, > a, is 
assumed (Fig. 9). Then, in agreement with the results of Sec. 3, we obtain: 
Zone I. The same solution is obtained as in Zone I of the first example. It is 
given by the Eqs. (4.2)-(4.5). 
Zone II. The equation of motion (3.1) has the form 
t 
As Py i ae 


(4.28) dv» Vo = 2 
1+-a,(t—%) : @ l+a,(t—1) 


dt 


-- 
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with the initial condition (3.2) 


(4.29) CA air a = A . 


On integrating, we obtain: 
t 
{1} 
v(t) = ° 


Q [-La,¢—t)) 
&=&, for t<t,.. (q, is given by the Eq. (4.5)), 


(4.31) | E> = Ee i for ii Stain 
fe 


(4.30) 


ay 

(4.32) ; 02(x, t) = Oo,(x)— Ao (x, t), 
where 

e P 0625 Ov,\ 
(4.33) Ao2(x, t) an i (Gee Fe 

1+a,(t—T) 
hence, after performing the calculations necessary: 

is 
(4.34) | ( ~3] 

(x, t) = Po| a2 eagtior) 
1 2(¢—2)[/-+a,(t—t,)]—(2t7—D)tay 
as pe (/+-a.(t—t,) 2 [x—/—a,(t—1))]}, 
a A ae 
Soa (! x] ; 

(4.35) i 

€ 

0, =—P, “(1 Hu for’) pee” aa, 
The stress jump is: 
(4.36) gre Poll ) eee 
2t Ay 


If ag = 0, we have oc, = 0, and the end is free. (The case of a, = 00 was discussed 
in the foregoing example. Now, in view of a, > a, this limit case cannot be inve- 
stigated). If a, = a,, we have o, = 0». 

Zone III, From the Eq. (4.28), for P(t) = 0 and the initial condition (3.4), we 
obtain: 
EAs 1 
4.37 AG ae ¢ 
( ) v3( ) 20 I+-a,(t—t,) 
€3 1S expressed by equations analogous to (4.31) 


(4.38) 03(X, t) = a2Po 


j 
\ 2[/+-a2(t—7)] 25) lt--ae(t=%) F 
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5. Conclusion 


The above solutions are obtained, as in [2], in a closed form owing to the as- 
sumption of rigid unloading characteristic, which leads, as indicated above, 
to certain qualitative changes in the process, but presents a relatively good quanti- 
tative approach to the reality. The combination of a great variety of cases with 
space restrictions placed on the present paper necessitates on discussing, only the 
most important cases. The rest may be solved by using the same methods. The 
assumption of a modified rigid unloading characteristic makes it possible to ob- 
tain the solution also, for the three dimensional states of stress for spherical and 
cylindrical waves, which is very simple in relation to the known numerical 
solutions. This will be the object of separate paper. 
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Sitwesze7/e nile 


FALA ODCIAZENIA DLA CIALA O SZTYWNEJ CHARAKTERYSTYCE ODCIAZENIA 
W OSRODKU WARSTWOWYM 


W pracy podano Scisle i zamkniete rozwiazania dla rozprzestrzeniania sie pla- 
styeznej fali odciazenia w oSrodku warstwowym dla ciata o sztywnej charaktery- 
styce odciazenia, przy granicy sprezystosci r6wnej zeru lub rdéznej od zera. Gdy 
o, #0, to zaklada sie réwniez sztywne obciazenie do granicy sprezystosci. 

W przypadku gdy a, = 0 otrzymuje sie rozwiazania otrzymane w pracy [2] dla 
odbicia sie fali od sztywnej Sciany przy o, = 0. 

W przypadku osrodka warstwowego io, 4 0 zjawisko rozprzestrzeniania si¢ fal 
ulega komplikacji. Ze wzgledu na liczne mozliwe warianty przypadkéw w zalez- 
nosci od przyjetych parametréw przedyskutowano zasadnicze sposréd nich. 

Rozwiazania dla pozostatych przypadkéw buduje sie analogicznie. Sztywna 
charakterystyka odciazenia daje pewne jakosciowe zmiany w obrazie rozprzestrze- 
nienia sie fal, wyniki iloSciowe otrzymuje sie jednakze bliskie rzeczywistosci. Poza 
tym odpada blad zwiazany z przyblizeniami matematycznymi, gdyz rozwiazania 
otrzymuje si¢ w postaci zamknietej. 
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Peswme 


BOJIHA PASIPY3KM JIA TEJIA, OBJIAJAIOWETO XKECTKOM 
XAPAKTEPMCTUKOM PASIPY3KH B CJIOMCTOM CPEJIE 


Jlaerca TOUHOe HM 33MKHYTOeC pelIeHMe AIA pacmpocrpaHeHuA WiacTMuecKoH 
BOJIHBI pasrpy3Ku B CNOMCTON cpeye WA Tema, oOMayaroujero »KECTKOM XapakTepu- 
CTHKOM pasrpy3KH, pH mpeyeme yupyroctw paBHOMYy HYyJIKO WIM He paBHOMY 
nyo. Korga o, #0 mpefmonaractca orcyrcrBue yupyrux jedopmaynit. 

B cmyuae ecu a,—=oo, MoyuaeTcA pelleHve, NouyueHHoe ye B padoTe 
[2] 1H oTpaxKeHHA BOJHBI OT 2KCCTKOM CTeHLI pH o, = 0. 

B cyyuae cyoucrol cpeqbi u o, #0 pacmpocrpaHeHue BONHbI yCIOXKHACTCA. 
Iipuuumaa BO BHHM@HHe MHOPOUNCJICHHbIC, BO3MO2%K BIC BapHaHTbI culyuaeB, B 
3aBHCUMOCTH OT IIPHHATBIX NapamMeTPOB, MpOBORMTCA WHCKycCaA MPHHIWWHasb- 
HbIX BapHaHTOB. 

PemieHue JIA OCTaIbHbIX CilyuaeB COCTaBIIAeTCA aHaworMuHoO. )KecTKad xapak- 
TepHCTHKA Pasrpy3KH JjacT HCKOTOPbIC KAUeECTbCHHbIe H3MCHCHHA B KapTHHe pac- 
MpOcrpaHeHHA BOJIH, OJHAKO KOJIMUCCTBCHHBIe Pe3yJIbTaTHI OM3KU WelCTBUTeIIb- 
HOocTH. Kpome Toro, MOXKHO H30@KAaTb OUWIMOKU, CBH3aHHO C MaTemMaTHUeCKUMU 
MpHOMWKEHHAMM TaK KaK PeWICHHA MOyuaIOTCA B 3aMKHYTOM BHI. 
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